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Nonlinear  dynamical  models  arise  from  both  theoretical  research  and  real  applications, 
including  robotics,  aircraft  control,  adaptive  control,  chemical  process,  numerical  algorithm, 
and  power  systems.  The  most  important  analysis  for  nonlinear  systems  is  the  determina- 
tion of  stability.  Only  critical  nonlinear  systems  need  further  investigation,  and  up  to  now, 
stability  analysis  has  been  predominantly  performed  using  center  manifold  and  bifurcation 
theories  in  combination  with  Lyapunov  techniques.  When  a  dynamic  system  is  locally 
asymptotically  stable  (LAS  in  short),  a  Lyapunov  function  must  exist.  Many  system  prop- 
erties, such  as  stability  region,  convergence  rate,  and  robustness,  are  associated  with  and 
can  be  estimated  using  Lyapunov  functions.  The  work  on  the  determination  of  the  local 
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stability  of  critical  nonlinear  systems,  or  the  construction  of  Lyapunov  functions  when  they 
are  LAS,  is  far  from  complete. 

In  this  dissertation,  critical  nonlinear  systems  are  studied  using  normal  form  theory. 
Once  a  critical  nonlinear  system  is  transformed  into  a  normal  form,  the  dynamics  of  the 
critical  states  are  decoupled  locally  from  those  of  the  stable  states.  Similar  to  results  derived 
using  center  manifold  theory,  the  stability  of  the  entire  nonlinear  system  is  determined  by 
that  of  critical  dynamics.  The  relationship  between  normal  form  and  center  manifold 
theories  is  then  explored.  If  the  critical  dynamics  is  LAS,  and  a  Lyapunov  function  is 
available,  a  Lyapunov  function  for  the  entire  system  can  be  easily  constructed. 

The  convergence  is  not  always  guaranteed  for  normal  forms.  Construction  of  Lya- 
punov functions  for  purely  critical  nonlinear  systems  is  only  known  in  a  few  special  cases. 
Therefore,  critical  nonlinear  systems  in  two  special  cases  are  studied  intensively  using  ap- 
proximate normal  forms.  The  explicit  stability  criteria,  as  well  as  explicit  algorithms  for 
construction  of  Lyapunov  functions,  are  developed.  Methods  for  simplification  of  Lyapunov 
functions  are  also  proposed.  The  construction  of  Lyapunov  functions  requires  deriving  (ap- 
proximate) normal  forms  and  the  associated  transformations,  and  is  very  computationally 
intensive.  The  issue  of  reducing  the  complexity  of  computing  normal  forms  and  the  associ- 
ated Lyapunov  functions  is  also  studied.  Critical  nonlinear  systems  with  multiple  critical 
eigenvalues  are  examined,  and  some  sufficient  conditions  for  determining  stability  and  con- 
structing Lyapunov  functions  are  obtained. 

This  dissertation  deals  with  continuous-time  systems.  All  of  the  results  are  applicable 
to  discrete- time  systems  with  slight  modifications. 
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CHAPTER  1 
INTRODUCTION 

Nonlinear  dynamical  models  arise  from  both  theoretical  research  and  real  applications, 
including  robotics,  aircraft  control,  adaptive  control,  chemical  processes,  numerical  algo- 
rithms, and  power  systems.  The  most  important  analysis  for  nonlinear  dynamical  systems 
is  the  determination  of  stability.  Linear  system  theory  has  be  well- developed  in  last  decades. 
To  a  certain  extent,  the  local  stability  of  a  nonlinear  system  can  be  completely  characterized 
by  its  linearization  at  the  equilibrium  point  of  interest.  The  local  stability  of  a  nonlinear 
system  cannot  be  determined  in  the  case  that  the  linearization  possesses  at  least  one  criti- 
cal eigenvalue  while  all  other  eigenvalues  are  stable.  This  dissertation  will  study  nonlinear 
dynamical  systems  whose  linearizations  exhibit  such  a  property.  The  local  stabilizability 
of  nonlinear  control  systems  possessing  linearly  uncontrollable  critical  eigenvalues  will  also 
be  discussed.  The  emphasis  will  be  on  the  construction  of  Lyapunov  functions  for  critical 
nonlinear  systems  which  are  known  to  be  locally  asymptotically  stable  (LAS).  This  chapter 
serves  as  an  introduction  to  the  dissertation.  In  Section  1.1,  we  present  much  of  the  back- 
ground material  which  will  be  used  throughout  this  dissertation.  Section  1.2  gives  a  brief 
literature  review.  The  organization,  along  with  the  major  contributions,  of  this  dissertation 
is  presented  in  Section  1.3. 

1.1    Problem  Formulation 

In  this  section,  we  present  much  of  the  background  material  that  will  be  used  throughout 
this  dissertation. 
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A  map  g  :  ft  C  Rm  -*  Rn  is  said  to  be  a  map  of  class  Ck  (or,  simply  Ck ),  denoted 
as  g  €  Cfc,  for  some  nonnegative  integer  k  if  its  partial  derivatives  of  up  to  order  A;  with 
respect  to  its  variables  exist  and  are  continuous  in  ft.  If  we  can  always  choose  k  —  +00, 
then  we  say  that  g  is  a  map  of  class  C°°  (or,  simply  C°°).  Furthermore,  g  is  said  to  be 
analytic  (denoted  as  C")  if  it  is  C°°  and  for  each  point  x0  €  ft  Taylor  series  expansion  of 
g  at  £0  converges  to  g(x)  for  all  a;  €  ft- 

The  goal  of  this  dissertation  is  to  study  the  local  asymptotic  stability  (LASY  in  short) 
of  nonlinear  dynamical  systems  of  form 

*=/(*),    xeRn  (1.1) 

where  /  :  ft  —  Rn  is  a  sufficiently  differentiate  (or  Ck  with  k  sufficiently  large)  map  from 
a  domain  ft  C  Rn  into  Rn . 

An  equilibrium  solution  (or  equilibrium  point,  or  simply  equilibrium)  of  (1.1)  is  a  point 

xo  €  ft  such  that 

/(x0)  =  0, 

i.e.,  a  solution  which  does  not  change  in  time. 

For  convenience,  we  state  all  the  definitions  and  theorems  for  the  case  when  the  equi- 
librium point  is  at  the  origin  of  Rn;  i.e.,  xQ  =  0.  There  is  no  loss  of  generality  in  doing 
this  since  any  equilibrium  point  can  be  shifted  to  the  origin  via  a  change  of  variables  [37]. 
Occasionally,  the  term  "zero  solution"  is  used,  instead,  to  signify  that  the  origin  is  the 
equilibrium  point  of  interest. 
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The  linearization  of  (1.1)  is  defined  as 

x=  Ax,  (1.2) 

where  A  =  Df(0)  is  the  first  order  derivative,  or  the  Jacobian  matrix  of  f(x)  :  Df(x)  = 
df(x)/dx,  evaluated  at  the  origin.  The  definition  of  eigenvalues  or  eigenvectors  is  well- 
known.  The  spectrum  of  an  n  x  n  matrix  A  is  the  set  of  eigenvalues  of  A,  and  denoted 
as  a(A)  =  {Ai,A2,---,An}  where  A,  needs  not  to  be  different  from  Xj.  Assume  A  is  an 
eigenvalue  of  the  n  x  n  matrix  A.  The  eigenvalue  A  is  said  to  be  hyperbolic  if  its  real  part  is 
nonzero  (denoted  as  Re(X)  ^  0),  stable  if  its  real  part  is  negative  (denoted  as  Re(X)  <  0), 
unstable  if  its  real  part  is  positive  (denoted  as  Re(X)  >  0),  and  critical  if  its  real  part  is 
zero  (denoted  as  Re(X)  =  0).  A  nonlinear  system  of  form  of  (1.1)  is  said  to  be  critical  if 
its  linearization  only  possesses  both  stable  and  critical  eigenvalues,  and  purely  critical  if 
its  linearization  only  possesses  critical  eigenvalues.  The  equilibrium  point  x0  =  0  is  said 
to  be  hyperbolic  if  none  of  the  eigenvalues  of  A  =  Df(0)  have  zero  real  part.  Similarly,  a 
nonlinear  control  system 

x  =  f(x,u) 

is  said  to  be  critical  if  the  linearization  about  x  =  0;  u  =  0  has  at  least  one  critical  eigenvalue 
that  is  not  controllable  while  all  other  eigenvalues  are  stabilizable. 

The  term  "stability"  has  different  definitions  in  different  areas.  In  this  dissertation, 
we  study  the  stability  of  nonlinear  systems  in  the  sense  of  Lyapunov.  Let  x(t)  denote  the 
solution  (trajectory)  of  the  system  (1.1)  which  at  the  initial  time  i0  =  0  passes  through  the 
initial  point  Xo  =  x(0). 

Definition  1  ([37],  [28])  The  equilibrium  point  x  =  0  of  (1.1)  is 
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•  stable  (in  the  sense  of  Lyapunov),  if  for  each  £  >  0,  there  exists  a  6(e)  >  0  such  that 

||as(0)||<«=>||x(OH<e,V«>0. 

•  unstable  (in  the  sense  of  Lyapunov),  if  it  is  not  stable. 

•  asymptotically  stable  (in  the  sense  of  Lyapunov),  if  it  is  stable  and  8  can  be  chosen 
such  that 

||i(0)||<£=>  lim  x(t)  =  Q. 

t — ►  OO 

00<> 

This  definition  is  mathematically  very  tidy.  However,  this  definition  does  not  provide 
us  with  a  method  for  determining  whether  or  not  a  given  equilibrium  point  is  stable.  The 
method  of  Lyapunov  can  often  be  used  to  determine  the  stability  of  equilibrium  points. 

Theorem  1  ([37],  [28])Consider  the  nonlinear  system  (1.1).  Let  x  —  0  be  an  equilibrium 
point  of  (1.1)  and  let  V  :  D  C  ft  R  be  a  C1  function  defined  on  some  neighborhood  D  of 
x  —  0  such  that 

V(0)   =    0     and     V(x)  >  0  Vj  ^  0 
V(x)    <    0     in  £-{0}. 

Then  the  equilibrium  point  x  =  0  is  locally  stable.  Moreover,  if 

V(x)  <  0    in    D  -  {0} 


then  the  equilibrium  point  x  =  0  is  locally  asymptotically  stable  (LAS  in  short). 


OOO 
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We  refer  to  V  as  a  Lyapunov  function.  When  the  origin  x  =  0  is  LAS,  we  are  often  in- 
terested in  determining  how  far  from  the  origin  can  the  trajectory  x(t)  be  and  still  converge 
to  the  origin  as  t  approaches  oo.  This  gives  rise  to  the  definition  of  the  region  of  attraction 
(stability  region,  region  of  asymptotic  stability,  domain  of  attraction,  or  basin).  The  region 
of  attraction  is  defined  as  the  set  of  all  points  x0  such  that  limt_>0o  x(t)  =  0  if  x(0)  =  xq. 
Finding  the  exact  region  of  attraction  analytically  might  be  difficult,  or  even  impossible. 
However,  Lyapunov  functions  can  be  used  to  estimate  the  region  of  attraction.  A  nonlinear 
system  is  said  to  be  globally  asymptotically  stable  (GAS  in  short)  if  and  only  if  its  region 
of  attraction  is  the  entire  space.  The  following  theorem  can  be  used  to  determine  whether 
or  not  a  nonlinear  system  is  GAS. 

Theorem  2  (Barbasin-Krasovskii  Theorem)  ([37],  [28])  Let  the  nonlinear  system  (1.1)  is 
defined  over  the  entire  space  Rn.  Let  V  :  Rn  ->  R  be  a  C1  function  such  that 

V(0)  =  0    and    V(x)  >  0  Vx  #  0 
||x||  — »■  oo  =>■  V(x)  —►  oo 
V(0)  =  0    and    V(x)  <  0  Vx  ^  0 

then  the  equilibrium  point  x  =  0  is  GAS.  OOO 

Remark  1  As  mentioned  before,  equilibrium  solutions  are  solutions  which  do  not  change  in 
time.  A  necessary  condition  for  x  =  0  to  be  LAS  in  Q  (or  GAS  in  Rn)  is  that  x  =  0  is  the 
only  equilibrium  point  in  Cl  (or  in  Rn). 

Example  1  (Lyapunov  function) 
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Consider  the  following  system: 

x    =    y  —  x3 

y    =  -x-x2y-y3. 

It  is  easy  to  verify  that  this  system  has  only  one  equilibrium  point  at  (x,y)  =  (0,0).  Let 
V(x,  y)  =  \{x2  +  y2).  Clearly  V(0, 0)  =  0,  V(x,  y)  —  oo  as  x2  +  y2  -  oo,  and  V(x,  y)  >  0 
V(x,y)  ^  0.  Then 

V  =  -x2y2  -  x4  -  y4  <  0  V(i,y)^0. 

Therefore  the  origin  is  globally  asymptotically  stable.  OOO 
The  linear  system  theory  has  been  well  developed  over  last  decades.  To  a  certain  extent, 
the  stability  and  the  instability  of  a  nonlinear  system  can  be  completely  characterized  by 
its  linearization. 

Theorem  3  [SOjConsider  the  system  (1.1).  Let  A  =  Df(0).  If  all  of  the  eigenvalues  of  A 
have  negative  real  parts,  then  the  origin  is  an  LAS  equilibrium  point.  If  A  possesses  any 
eigenvalue  with  positive  real  part,  then  the  origin  is  unstable.  If  A  only  possesses  both  stable 
and  critical  eigenvalues,  the  stability  of  (1.1)  cannot  be  deduced  from  the  linearization. 

For  later  purposes,  we  need  the  following  two  theorems  which  give  criteria  for  deter- 
mining local  stability  of  nonlinear  systems  in  two  special  cases. 

Theorem  A  [49]  A  one  dimensional  analytic  system  of  form 


x  =  aix  +  a^x2  +  a^x3  + 
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is  LAS  if  and  only  if  a,k  <  0  and  k  is  odd  where  k  =  min{t  |  a,  ^  0}.  i.e.,  the  first  nonzero 
number  in  the  sequence  {a\,  a2, a3,  •  •  • ,  }  is  negative  and  the  subscript  is  odd.  000 

This  theorem  can  be  easily  proved  using  the  Lyapunov  function  V(x)  =  \x2. 
Now  consider  a  two-dimensional  analytic  system  of  form  with  u  €  R+ 

x    =    -uy  +  a2ox2  +  anxy  +  a02y2 

+a30x3  +  a2lx2y  +  a12xy2  +  aQ3y3  +  0(\x,  y\4) 

y    =    ux  +  b20x2  +  buxy  +  b02y2 

+630X3  +  b2lx2y  +  b12xy2  +  b03y3  +  0{\x,  y\4). 

After  a  proper  change  of  coordinates  (x,  y)  -*  (f ,  77),  we  obtain 

£   =    -un  +  f^ie  +  ^m^-n^  +  Oi^nl2^2)  (1.5) 
i=i 

V    =   ^  +  E(^  +  r?2),(m2l+1ry  +  n2l+10  +  O(|^J?r+2).  (1.6) 
!=i 

Theorem  5  ([6],  [27])The  system  (1.5),  (1.6)(or  equivalently  (1.3),  (1-4)  )  is  LAS  if  and 

only  if  m2k+i  <  0  where  k  =  min{i  |  m2i+i  ^  0},  i.e.,  the  first  nonzero  number  in  the 

sequence  {m3,ms,  m.7,  •  •  ■ ,  }  is  negative.  000 

This  theorem  can  be  easily  proved  using  the  Lyapunov  function  V(£,  n)  =  \{£2  +  n2)  for 
the  system  (1.5),  (1.6).  As  we  will  see  later,  the  system  (1.5)  and  (1.6)  is  one  of  the  normal 
forms  of  (1.3)  and  (1.4).  For  a  given  system,  the  computation  of  m2i+i  is  quite  involved 
but  systematic.  For  example, 


(1.3) 


(1-4) 
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+  tt—  [<*n(a2o  +  002)  -  hi{ho  +  &02)  -  G20&20  +  ^02^02]  • 
loo; 

1.2    Brief  Literature  Review 

Stability  (or  stabilizability)  of  critical  nonlinear  (control)  systems  is  a  classical  problem 
in  both  theoretical  research  ([6],  [38],  [53],  [16],  [15])  and  variety  of  applications  ([41], 
[43],  [40],  [46],  [52],  [45]).  Many  different  approaches  are  used  to  study  nonlinear  stability 
(or  stabilizability)  analysis.  In  particular,  center  manifold  theory  is  used  in  [3],  [11],  [6], 
[39],  and  [7].  In  [1],  [2],  [25],  [32],  [33],  [34],  [35],  [41],  [43],  and  [52],  bifurcation  theory 
is  employed.  Even  though  center  manifold  and  bifurcation  theories  are  successfully  used 
in  stability  (or  stabilizability)  analysis;  these  techniques  are  found  (up  to  the  date  this 
dissertation  is  finished,  and  to  author's  best  knowledge)  not  very  helpful  for  the  construction 
of  Lyapunov  functions  when  a  critical  nonlinear  system  is  LAS. 

In  this  dissertation,  we  propose  to  study  critical  nonlinear  systems  using  Poincare  nor- 
mal form  theory.  The  normal  form  method  is  a  technique  of  change  of  coordinates,  and  is 
capable  of  bringing  out  the  essential  features  of  nonlinear  systems.  The  technique  of  using 
change  of  coordinates  in  nonlinear  system  analysis  is  not  new.  Lyapunov  already  used  tech- 
nique  of  change  of  coordinates  in  stability  analysis  of  nonlinear  systems  in  his  landmark 
Ph.D.  dissertation  [42]  in  1892.  Different  normal  form  techniques  are  discussed  in  much 
of  literature,  such  as  [4],  [5],  [8],  [9],  [14],  [19],  [20],  [17],  [18],  [21],  [22],  [36],  [31],  [50]. 
Roughly  speaking,  the  linearization  of  affine  nonlinear  systems  [29],  [47]  also  belongs  to  the 
area  of  normal  form  techniques.  In  [6],  and  [38],  the  Poincare  normal  form  is  employed 
for  the  stability  analysis  of  critical  nonlinear  systems.  Behtash  and  Sastry  [7]  examined 
nonlinear  (control)  systems  using  both  normal  form  and  center  manifold  theories.  However, 
the  power  of  normal  form  techniques  has  not  yet  been  fully  recognized  and  exploited.  The 
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search  to  discover  other  potential  merits  of  normal  form  techniques  in  study  of  nonlinear 
(control)  systems  motivates  much  of  the  work  presented  in  this  dissertation. 

When  a  dynamical  system  is  LAS,  the  existence  of  Lyapunov  functions  is  guaranteed. 
For  linear  systems  or  nonlinear  systems  whose  linearization  possesses  no  critical  eigenvalues, 
the  algorithm  for  the  construction  of  Lyapunov  functions  is  well  established.  Many  system 
properties  such  as  stability  region,  convergence  rate,  and  robustness  can  be  estimated  using 
Lyapunov  functions  (see  [26],  [12],  [13]  for  example).  Unfortunately,  there  is  no  general 
procedure  for  the  construction  of  Lyapunov  functions  for  critical  nonlinear  systems.  In 
[44],  Mees  and  Chua  constructed  a  Lyapunov  function  for  a  simple  2-dimensional  nonlinear 
system  with  a  pair  of  purely  imaginary  eigenvalues.  The  method  they  used  cannot  be 
generalized.  Fu  and  Abed  [24]  proposed  a  procedure  for  the  construction  of  Lyapunov 
functions  for  nonlinear  systems  in  two  special  cases.  The  procedure  is  based  on  a  lemma 
about  the  local  definiteness  of  a  class  of  bivariate  functions,  and  is  only  applicable  when  the 
local  stability  of  a  nonlinear  system  can  be  determined  using  only  the  first  through  the  third 
order  terms  without  necessity  of  the  consideration  of  the  higher  order  terms.  However,  with 
the  help  of  Poincare  normal  form  techniques,  the  algorithm  proposed  in  this  dissertation  is 
able  to  construct  Lyapunov  functions  for  nonlinear  systems  in  two  special  cases  no  matter 
how  many  terms  in  the  dynamics  are  necessary  to  determine  the  local  stability.  Other 
nonlinear  systems  with  more  than  one  (or  one  pair  of  )  critical  eigenvalue(s)  is  also  studied 
using  normal  form  techniques. 
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1.3    Outline  of  the  Dissertation 

Aside  from  the  first  and  last  Chapters,  this  dissertation  can  be  divided  into  two  parts. 
The  mathematical  preliminaries  are  introduced  in  Chapters  2  and  3,  while  most  of  results 
from  this  research  are  presented  in  Chapters  4  through  8. 

In  Chapter  2,  center  manifold  theory  is  introduced.  Center  manifold  theory  is  a  pre- 
vailing model  reduction  method  in  analysis  of  critical  nonlinear  systems.  We  include  the 
definition  of  center  manifold  and  three  fundamental  theorems  of  the  center  manifold  theory 
in  this  chapter.  Several  examples  are  given  to  illustrate  use  of  the  center  manifold  theory. 

Chapter  3  presents  the  introduction  to  normal  form  theory.  In  essence,  normal  form 
theory  is  a  change  of  coordinates  technique.  We  first  illustrate  the  concept  and  structure 
of  normal  forms  using  near  identity  of  changes  of  coordinates.  An  efficient  algorithm  for 
the  computation  of  normal  forms  follows.  By  introducing  and  studying  a  linear  operator 
in  the  space  of  homogeneous  maps  with  a  properly  defined  inner  product,  the  structure  of 
normal  forms  can  be  predicted.  The  space  of  homogeneous  maps  is  a  linear  one.  It  is  also 
possible  to  study  normal  forms  in  other  linear  spaces  once  a  correspondence  is  established 
between  the  two  linear  spaces. 

In  Chapter  4,  we  investigate  the  relationship  between  center  manifold  and  normal  form 
theories.  Center  manifolds  for  a  given  critical  nonlinear  system  need  not  be  unique.  How- 
ever, it  is  shown  that  all  center  manifolds  are  equivalent  in  some  sense.  When  normal 
form  theory  is  applied  to  a  critical  nonlinear  system,  the  system  in  normal  forms  exhibits 
a  decoupled-like  structure.  The  stability  of  the  entire  system  is  determined  by  that  of  the 
critical  dynamics.  This  conclusion  is  identical  to  the  results  derived  using  center  mani- 
fold theory.  The  relationship  between  center  manifold  and  normal  form  theories  is  then 
explored. 
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The  merits  of  normal  forms  are  exploited  further  in  Chapter  5.  We  develop  procedures 
for  construction  of  Lyapunov  functions  for  critical  nonlinear  systems  when  they  are  locally 
asymptotically  stable.  First,  a  general  procedure  is  provided.  Since  the  critical  dynam- 
ics are  independent  of  stable  dynamics  for  systems  in  normal  forms  (when  they  exist), 
Lyapunov  functions  for  the  entire  system  can  be  derived  by  just  summing  the  individual 
Lyapunov  functions  corresponding  to  critical  and  stable  dynamics,  respectively.  However, 
the  normal  form  procedure  is  not  always  convergent,  and  Lyapunov  functions  for  purely 
critical  systems  are  not  always  available.  Nevertheless,  we  are  able  to  develop  an  algorithm 
for  the  construction  of  Lyapunov  functions  for  critical  nonlinear  systems  in  two  specific 
cases  using  approximate  normal  forms. 

Compared  to  the  Lyapunov  functions  constructed  by  other  researchers  for  some  simple 
examples,  the  Lyapunov  functions  constructed  in  Chapter  5  are  quite  complicated  in  some 
sense.  In  Chapter  6,  an  algorithm  for  the  simplification  of  Lyapunov  functions  is  proposed. 
Several  examples  are  given  to  highlight  the  simplicity  of  these  refined  Lyapunov  functions. 

An  alternative  algorithm  (compared  to  the  one  given  in  Chapter  5)  is  proposed  in 
Chapter  7.  To  use  this  algorithm,  a  modified  (partial)  approximate  normal  form  is  sufficient 
for  the  construction  of  Lyapunov  functions.  This  will  significantly  reduce  the  complexity 
of  computing  normal  forms  and  the  associated  transformations.  The  simplification  of  the 
associated  Lyapunov  functions  is  also  addressed. 

Chapter  8  includes  several  independent  ideas.  We  apply  normal  form  techniques  to 
a  broader  class  of  critical  nonlinear  systems  than  those  considered  in  previous  chapters. 
Under  certain  hypotheses,  the  stability  of  such  systems  can  be  easily  determined  and  the 
associated  Lyapunov  functions  can  be  easily  constructed.  The  bulk  of  the  research  presented 
in  this  dissertation  only  treats  systems  in  the  continuous-time  domain.  Actually,  normal 
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form  techniques  are  applicable  to  discrete-time  systems  using  a  slight  modification.  An  brief 
introduction  to  use  of  normal  form  theory  in  the  discrete-time  domain  concludes  Chapter 
8. 

Finally,  in  Chapter  9,  we  summarize  the  dissertation  and  give  some  concluding  remarks. 


CHAPTER  2 

INTRODUCTION  TO  CENTER  MANIFOLD  THEORY 

The  center  manifold  theorem  in  finite  dimensions  can  be  traced  to  the  work  of  Pliss 
[48].  Other  valuable  references  include  [11],  [53],  [27],  and  [15]. 

Center  manifold  theory  provides  a  model  reduction  technique  for  analyzing  critical  non- 
linear dynamic  systems.  As  we  have  mentioned  before,  a  necessary  condition  for  a  nonlinear 
system  to  be  LAS  is  that  the  linearization  possesses  no  eigenvalues  with  positive  real  parts. 
A  nonlinear  dynamic  system  is  LAS  if  all  eigenvalues  of  the  linearization  have  negative  real 
parts.  The  only  case  when  the  local  stability  of  a  nonlinear  system  cannot  be  deduced 
from  the  linearization  is  when  the  linearization  possesses  both  stable  and  critical  eigenval- 
ues. In  this  situation,  nonlinear  terms  must  be  taken  into  consideration  for  determining 
the  stability  of  a  nonlinear  dynamical  system.  Center  manifold  theory  allows  us  to  deduce 
the  stability  of  a  nonlinear  system  by  studying  a  purely  critical  nonlinear  system  which  in 
general  has  a  reduced  dimension. 

The  definition  of  the  center  manifold  is  given  below  first. 

Definition  2  [53,  11]  (Invariant  Manifold)  A  set  S  €  Rn  is  said  to  be  a  local  invariant 
manifold  for  the  dynamical  system 

x=f(x),  xeRn  (2.1) 
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if  for  xQ  €  S,  the  solution  x(t)  of  (2.1)  with  x(0)  =  x0  is  in  S  for  \t\  <  T  where  T  >  0.  // 
we  can  always  choose  T  =  +oo,  then  we  say  that  S  is  an  invariant  manifold.  OOO 

Definition  3  [53,  11]  (Center  Manifold)  An  invariant  manifold  S  is  called  a  center  manifold 
for  (2.2)-(2.3)  if  it  can  be  locally  represented  as  follows 

S  =  {{xc,xs)eRncxRn'\x3  =  h{xc);   \x\<e;   ft(0)  =  0;   Dh(0)  =  0} 

for  some  e  >  0.  <>oc> 

Center  manifold  theory  has  several  applications  to  dynamical  systems.  In  this  research, 
center  manifold  theory  is  employed  for  determining  the  stability  of  dynamical  systems  of 
form 

xc  =  Acxc  +  fc(xc,  xs)  (2.2) 
xs    =    As  xs  +  fs(xc,  x3),  (2.3) 

where  xc  €  RUc,  x„  G  Rn° ,  A°  and  A3  are  constant  matrices  such  that  all  of  the  eigenvalues 
of  A°  have  zero  real  parts  while  all  of  the  eigenvalues  of  A3  have  negative  real  parts.  The 
vector  valued  functions  /c,  /,  are  at  least  Ck  (k  >  2)  with  /c(0,0)  =  0,  /,(0,0)  =  0, 

£/c(0,0)  =  0,  £>/.(0,0)  =  0. 

The  three  fundamental  theorems  of  the  center  manifold  theory  are  summarized  as  fol- 
lows: 

Theorem  6  [11](  Existence  of  Center  Manifolds)  There  exists  a  Ck  center  manifold  y  -  h(x) 
for  (2.2)-(2.3).    The  dynamics  of  (2.2)-(2.3)  restricted  to  the  center  manifold  is,  for  £ 
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sufficiently  small,  given  by  the  following  nc-dimensional  system 


Z=Act  +  fc(Z,  h(0),URnc-  (2-4) 


ooo 


Theorem  7  [11] (Reduction  Principle)  (I)  Suppose  the  zero  solution  of  (2.4)  is  stable  (or 
asymptotically  stable,  or  unstable).  Then  the  zero  solution  of  (2.2)-(2.3)  is  also  stable  (or 
asymptotically  stable,  or  unstable).  (II)  Suppose  the  zero  solution  of  (2.4)  is  stable.  Let 
(xc(t),  x3(t))  be  a  solution  of  (2.2)-(2.3)  with  (xc(0),  x,(Q))  sufficiently  small.  Then  there 
is  a  solution  of  (2-4)  such  that  as  t  — ►  +oc 

xc(t)  =  f(t)  +  0(e^*) 
x,(t)  =  *(«*)) +  0(«"7*) 

where  7  >  0  is  a  constant.  OOO 
Substituting  the  center  manifold  x,(t)  =  h(xc(t))  into  (2.3),  we  obtain 

Dh(xe)[Aexe  +  fc(xc,  h(xe))]  =  Aah{xc)  +  fs(xc,  h(xc)).  (2.5) 

The  center  manifold  y  =  h(x)  is  obtained  by  solving  equation  (2.5)  together  with  the 
conditions  h(0)  =  0,  Dh(0)  —  0.  A  closed  form  solution  to  the  PDE  (2.5)  is  often  not  easily 
found.  The  next  theorem,  however,  shows  that  the  center  manifold  can  be  approximated, 
in  principle,  to  any  degree  of  accuracy. 
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Theorem  8  [11] (Approximation  of  the  Center  Manifold)  Let  <f>  :  D  C  Rn<  —  R71'  be  a  C1 
mapping  with  <f>(0)  =  0  and  D<f>(0)  =  0.  Suppose  that  as  xc  —  0 


(M<j>)(xc)  =  D<p(xc)[Acxc  +  fc(xc,4>(xc))]  -  A"<f>(xc)  -  fs{xc,<t>{xc))  =  0(\\xe\\") 


for  some  q>l.  Then  as  xc  ->  0,  \\h(xc)  -  <i>{xc)\\  =  0(| |arc| I'').  OOO 

Remark  2  [11]  Center  manifolds  need  not  to  be  unique.  However,  there  are  some  points 
which  must  always  lie  on  the  any  center  manifold.  For  example,  any  center  manifold  must 
include  the  equilibrium  point.  OOO 

Remark  3  [11]  If  fc  and  fs  are  Ck,  (2  <  k  <  oc),  then  h  is  also  Ck .  If  fc  and  f3  are 
analytic,  then  in  general  (2.2)-(2.3)  does  not  necessarily  have  an  analytic  center  manifold, 
not  even  a  C°°  one.  OOO 

A  few  examples  are  provided  to  illustrate  the  use  of  center  manifold  theory. 
Example  2  (Derivation  of  the  center  manifold) 
Consider  the  system 

xc  -  xcxs  +  ax3c  +  bxcx]  (2.6) 
xs    =    -xs  +  cx\  +  dx\xs  +  ex\,  (2.7) 

where  a,  b,  c,  d,  and  e  are  real  constants. 

By  Theorem  6,  the  system  (2.6)-(2.7)  has  a  center  manifold  x3  =  h(xc).  Assuming 
h(xc)  =  h2x2c  +  h3x3  +  h4x*  +  h5x^.  +  •  •  •,  we  obtain 


(Mh)(xc)    =    (-c  +  h2)x2c  +  (h3-e)x3c 


+  (2h2  (h2  +  a)  +  h4-dh2)x4c 

+  (h5  -dh3  +  2  h2h3  +  3  h3  (h2  +  a))  x5c+0  (\\xc\\6)  . 

Hence,  (Mh)(xc)  =  0  (\\xcf)  if  we  take 

h(xc)  =  ex]  +  ex3c  +  (dc  -  2c(c  +  a))x4c  +  (de  -  2ce  -  3e(c  +  a))x5c  +  ■■■. 

By  Theorem  7,  the  dynamics  which  determines  the  stability  of  the  zero  solution  of  the 
system  (2.6)-(2.7)  is 


(c  +  a)?  +      +  (dc  -  2c(c  +  a)  +  be2)  £5  +  0 


Example  3  [11]  (Non-uniqueness  of  the  center  manifold) 


Consider  the  system 


OOO 


xr  =  -x: 


(2.8) 
(2.9) 


It  is  easy  to  verify  that  for  any  real  number  Ci,c2, 


xs  =  h(xc,ci,c2)  =  < 


c1exp(-^xc  2),    xc  >  0 


0  x,  =  0 


c2exp(-i;rc2),    xc  <  0 


satisfies  (2.5),  hence  is  a  center  manifold  for  the  system  (2.8)-(2.9). 
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Example  4  [51]  (Non-analyticity  of  the  center  manifold) 
Consider  the  system 

xe   =    -x\  (2.10) 
=    -xs  +  x2c.  (2.11) 

Suppose  that  the  system  (2.10)-(2.11)  has  a  center  manifold  which  is  analytic  at  the 
origin.  Then,  by  Taylor's  Theorem, 

oo 

h(xe)  =  £  hnxnc  (2.12) 

n=2 

in  some  neighborhood  of  the  origin.  The  coefficients  {hn}  are  determined  by  substituting 
(2.12)  into  (2.5).  We  have, 

h2  =  1 

^2n+i  =  0,    for    n  =  1,  2,  3,  •  •  • 
/in+2  —  nhn,    for    n  —  2,  4,  6,  •  •  •  . 

A  series  of  the  form  (2.12)  with  these  coefficients  is  not  convergent.  Therefore  the  system 
(2.10)-(2.11)  cannot  have  an  analytic  center  manifold.  OOO 

Example  5  (Nonlinear  systems  depending  on  parameters) 

Consider  the  system 

x    =   sx-x3  +  xy  (2.13) 

V   =    -y  +  y2-x2  (2.14) 
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where  £  is  a  real  constant  parameter. 

It  is  evident  that  the  origin  (0,0)  is  an  isolated  equilibrium  point  for  the  system  (2.13)- 
(2.14).  The  eigenvalues  of  the  Jacobian  matrix  corresponding  to  this  equilibrium  point 
are  -1  and  e.  Center  manifold  techniques  cannot  be  applied  to  (2.13)-(2.14)  directly  since 
the  sign  off  is  unknown.  However,  we  can  re-write  (2.13)-(2.14)  in  an  equivalent  form  as 
follows 

i  =  ex  -  x3  +  xy  (2.15) 
y  =  -y  +  y2-x2  (2.16) 
i    =    0.  (2.17) 

When  considered  as  an  equation  on  R3  the  ex  term  in  (2.15)  is  nonlinear.  Thus  {-1,0,0} 
is  the  set  of  eigenvalues  for  the  Jacobian  matrix  corresponding  (2.15)-(2.17)  at  the  origin 
(0,0,0).  Now  center  manifold  theory  can  be  applied.  Assume  a  center  manifold  for  the 
system  (2.15)-(2.17)  is 

y  =  h(x,  e)  =  a2ox2  +  anxe  +  a02e2  +  ••-, 

we  obtain 

(Mh)(x,  e)  =  (1  +  a20)x2  +  anxe  +  a02e2  +  0{\x,  e\3). 

Thus  y  =  h(x,e)  =  —x2  +  ■  ■  ■  is  a  2nd  order  approximation  to  center  manifolds.  Since 
for  any  value  of  e,  (0,0,  £•)  is  always  an  equilibrium  point  for  the  system  (2.15)-(2.17)  (see 
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Remark  2),  we  conclude  that  h(Q,e)  =  0.  The  dynamics  along  the  center  manifold  is 

e  =  *  -  2£3+ om  -im3) 

e   =  0. 

It  is  easy  to  conclude  now  that  the  system  (2.15)-(2.17)  is  LAS  if  e  <  0,  and  unstable  if 
e  >  0. 

As  we  remarked  before,  center  manifolds  need  not  to  be  unique  for  a  given  nonlinear 
system.  However,  it  will  be  shown  in  Chapter  4  that  all  of  the  center  manifolds  for  a  given 
nonlinear  system  are  equivalent  in  some  sense. 


CHAPTER  3  , 
INTRODUCTION  TO  NORMAL  FORM  THEORY 

This  chapter  is  concerned  with  normal  form  theory.  The  method  of  normal  forms  can 
be  traced  to  the  Ph.D  thesis  of  Poincare,  is  further  refined  in  [50].  Other  good  references 
include  [4],  [17],  [17],  [15],  and  [10].  In  Section  3.1,  the  concept  as  well  as  structure  of 
normal  forms  are  introduced.  Derivation  of  normal  forms  is  very  computationally  intensive. 
Efficient  algorithms  are  presented  in  Section  3.2.  In  Section  3.3,  we  give  a  rigorous  definition 
of  the  normal  form,  and  provide  some  observations  about  the  linear  operator  LkA.  Several 
examples  are  provided  in  this  chapter  as  an  illustration  of  the  use  of  normal  form  theory. 
A  brief  summary  of  this  chapter  is  given  in  Section  3.4. 

3.1    Introduction  to  Normal  Forms 

The  concept  as  well  as  structure  of  normal  forms  can  be  better  illustrated  by  applying 
a  sequence  of  near  identity  changes  of  coordinates  to  a  nonlinear  system. 

To  simplify  the  description  of  normal  forms  for  nonlinear  dynamical  systems,  it  is  first 
necessary  to  introduce  some  basic  notation  about  homogeneous  polynomial  maps  of  a  given 
degree.  Other  related  notation  will  be  defined  sequentially  as  we  go  along.  All  of  this 
notation  is  defined  over  the  fields  of  real  numbers,  but  is  also  valid  over  the  field  of  complex 
numbers. 

For  a  positive  integer  k  >  0,  a  monomial  of  degree  k  over  Rn  is  an  expression 

=  with    |a|  =  £>  =  * 
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where  a  =  [ax,  a2,  •  •  ■ ,  an]  (with  integer  a,-  >  0)  is  a  multi-index. 

A  homogeneous  polynomial  function  of  degree  k,  f  :  Rn  -*  R  is  defined  as  the  linear 
combination  of  all  monomials  of  degree  k  over  Rn,  i.e., 

/=  L  6«*a 

where  ba  £  R.  A  vector-valued  polynomial  function  (map)  f  :  Rn  —>  Rm  is  said  to  be 
homogeneous  of  degree  fc  if  every  entry  is  a  homogeneous  polynomial  function  of  degree  k. 

In  this  dissertation,  a  homogeneous  polynomial  map  (function)  of  degree  k  will  be 
denoted  as  fk  with  the  superscript  indicating  the  degree.  Notice  that  zero  is  a  homogeneous 
polynomial  map  (function)  of  any  degree. 

The  space  of  homogeneous  polynomial  maps  from  Rn  to  Rm  of  degree  k  is  defined  as 

H(n,m,k)  =  {fk  :  Rn  -  Rm}. 

Therefore,  H(n,m,k)  is  a  linear  space  with  dim  H(n,  m,  k)  =  m(n+j^_1)  =  rn^Zly'  ■  An 
inner  product  (■,•)  on  the  space  H(n,m,k)  is  introduced  in  [4]. 

Let  A  be  an  n  x  n  real  constant  matrix.  For  each  k  >  2,  the  linear  operator  L\  : 
H(n,  n,  k)  -+  H(n,  n,  k)  is  defined  by 

(LkAPk)(y)  =  DPk{y)Ay-APk{y)  =  [Ay,  Pk{y)}. 

This  operator  will  be  used  in  constructing  normal  forms. 
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Let  R(n,n,k)  denote  the  range  of  LkA,  and  C(n,n,k)  any  complementary  (not  neces- 
sarily the  orthogonal )  subspace  of  R(n,  n,  k)  in  H(n,  n,  k).  This  sets  up  the  decomposition 

H(n,n,k)=  R(n,n,k)  ®  C(n,n,k).  (3.1) 

Now  we  are  in  the  position  to  introduce  normal  form  theory.  Without  loss  of  generality, 
consider  an  analytic  nonlinear  dynamical  system  expanded  as  a  series  of  homogeneous 
polynomial  maps  as  follows: 

i=  /(»)  =  /»(*)  +  f2(x)  +■■■  +  fk~\x)  +  fk(x)  +  fk+1(x)  +  ■■■,       x  e  R\  (3.2) 

where  /'(•)  G  H(n,n,i).  Here  f\x)  =  Ax  with  A  =  Df{0). 

For  Pk(-)  £  H(n,n,k),  consider  now  the  effect  of  a  near  identity  change  of  coordinates 
of  the  form 

x  =  y  +  Pk(y).  (3.3) 
Differentiating  both  sides  of  (3.3)  with  respect  to  time  and  using  (3.2)  yields 

(I+DPk(y))y=f(y  +  Pk(y)).  (3.4) 

Note  that  (/  +  DPk(y))  is  always  invertible  in  a  sufficiently  small  neighborhood  of  the 
origin,  and 

(/  +  DPk(y)Tx  =1-  DPk(y)  +  0(\y\2k-2),  (3.5) 
where  0(\y\2k~2)  denotes  combination  of  terms  of  degree  (2k  -  2)  or  higher. 
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Equations  (3.4),  (3.5)  and  (3.2)  yields  the  dynamics: 

i/=  g(y)  =  gl(y)  +  gl(y)  +  •  •  • + 9k~Hv)  +  /*(»)-  +  •  •  ■ ,  (3-6) 

where  yj(-)  =  /»'(•)  for  i  =  1,2,  •  •  -  1.  In  general,  for  i  >  k  +  1,  flj(-)  ^  /*"(■)  in  the 
equation  (3.2).  Here,  the  subscript  'it'  corresponds  to  the  fact  that  the  transformation 
Pk  £  H(n,  n,  k).  The  term  of  degree  k  in  (3.6)  is  given  by 

gkk(y)  =  fk(y)  -  [Ay,  Pk(y)]  =  fk(y)  -  (LkAPk)(y).  (3.7) 

This  suggests  introducing  the  linear  operator  LkA\  H(n,n,k)  -»  H(n,n,k)  for  each  A;  >  2. 

If  fk  €  R(n,n,k),  we  can  choose  an  appropriate  Pk(-)  such  that  (LkAPk){y)  =  fk{y), 
and  consequently  completely  cancel  out  all  terms  of  degree  k  in  (3.6).  In  the  more  gen- 
eral cases  where  fk  $  R(n,n,k),  we  can  always  pick  up  an  appropriate  Pk{-)  such  that 
fk(y)  _  [Ay,Pk(y)]  e  C{n,n,k).  The  terms  of  degree  k  which  cannot  be  removed  should 
be  considered  in  some  sense  as  the  essential  part  for  the  dynamical  behavior  of  (3.2)  near 
the  origin. 

In  order  to  implement  the  normal  form  transformation,  apply  a  sequence  of  changes 
of  coordinates  (3.6)  to  the  system  (3.2)  iteratively,  starting  with  k  =  2.  At  each  step, 
an  appropriate  Pk  is  chosen  such  that  g^iy)  -  [Ay,Pk{y)]  €  C(n,n,k).  After  a  finite 
number  of  iterations,  the  system  (3.2)  is  said  to  be  put  into  an  approximate  normal  form. 
Clearly,  the  process  of  "normalization"  is  intrinsically  local  in  the  sense  that  the  change  of 
coordinates  is  only  valid  in  some  neighborhood  of  the  origin.  The  following  normal  form 
theorem  is  obtained  from  the  above  discussion. 
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Theorem  9  (Normal  Form  Theorem) [4],  [6]  For  the  system  (3.2),  given  the  sequence  of 
decompositions  (3.1)  and  any  finite  positive  integer  k,  there  exists  a  series  of  changes  of 
coordinates,         =  x,  +  P'(x;),  X\  =  x,  xjt  =  y,    i  =  2,  3,  P'  £  H(n,n,i),  and 

a  neighborhood  ft  of  the  origin,  such  that  for  y  £  SI,  the  system  (3.2)  is  transformed  to  a 
system  of  the  form 

y  =  g\y )  +  g\v)  +  •  ■  •  +  9k(y)  +  0(\\y\\k+1),  (3.8) 

where  g\y)  =  Ay  and  g{(y)  £  C(n,  n,  i)  for  k  =  2, 3,  •  ■  ■ ,  k.  OOO 

Remark  i  Clearly  the  structure,  as  well  as  actual  expression,  of  the  normal  form  (3.8)  of 
(3.2)  depends  on  A.  Even  if  A  is  fixed,  the  normal  form  is  still  not  unique,  but  depends  on 
the  choice  of  complementary  subspaces  C(n,  n,  i),  as  well  as  other  factors.  A  more  rigorous 
definition  of  the  normal  form  will  be  given  later.  OOO 

Remark  5  As  long  as  the  number  of  the  iterations  of  the  changes  of  coordinates  is  finite, 
the  normal  form  (3.8)  of  (3.2)  must  be  convergent  in  some  neighborhood  Q.  of  the  origin. 
If  the  sequence  of  changes  of  coordinates  is  taken  to  infinity,  the  resulting  sequence  of 
transformations  may  not  be  convergent.  In  the  case  that  the  resulting  normal  form  process 
is  convergent,  and  f  €      ,  the  resulting  full  normal  form  is  given  by: 

y  =  g(y)  =  g\y)  +  g2(y)  +  93(y)  +  ■■■  (3-9) 


where  g1(y)  =  Ay  and  gk  €  C(n,  n,  k)  for  k  =  2, 3, 
Example  6  (Derivation  of  Normal  Forms) 


OOO 
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Consider  the  system 


x    =    xy  +  x" 


y   =    -y-2x  , 


(3.10) 
(3.11) 


where  A  = 


0  0 
0  -1 


For  terms  of  degree  two,  we  have 


it    .  \ 


#(2,2,2)  =  span  < 


IV  0  ) 


I  \ 

xy 


\  °  / 


/  \ 


V  0  / 


x2 


I  \ 


(      \  1 


\V   I  J 


>  . 


We  want  to  compute  L2A(H{2, 2, 2)).  We  do  this  by  applying  the  linear  operator  L\  to  each 
basis  element  on  #(2,2,2).  These  are: 


\  °  / 
/  \ 

xy 


\  °  / 


V  0  / 

/  \ 

o 


x2 

/  \ 
0 


\xy ) 

I  \ 

0 


I 


2x  0 
0  0 
y  x 
0  0 
0  2y 
0  0 
0  0 
2x  0 
0  0 
y  x 
0  0 
0  2y 


\  I       \  I 
0 


\ 


\  I       \  I 
0 


\-y ) 


\  I    \  t 

0 


\-y ) 


\  I    \  ( 

0 


I  \-y  J  \ 
\  (    \  t 

0 


\-y  J 

\  I    \  ( 

0 


\-y ) 


0  0 
0  -1 
0  0 
0  -1 
0  0 
0  -1 
0  0 
0  -1 
0  0 
0  -1 
0  0 
0  -1 


V  0  / 

\  (  \ 

xy 


\  °  / 
\  /  \ 

v2 


\ 0  / 


\  /  \ 

0 


x2 
\  ) 


0 


\xy ) 

\  I  \ 

0 


(  \ 

0 


v0/ 


=  -2 


/  \ 
xy 


\  °  / 


V  0  / 


/  \ 

o 


x2 
\  I 


I  \ 
0 


V0/ 


0 
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it  \ 


Therefore  il(2,2,2)  =  span  < 


IV  u  J 


(  \ 


v°  / 


/  \ 


\x ) 


I  w 

o 


>  ,  C(2, 2, 2)  can  be 


V 

/  \ 

< 

0 

V 0  / 

[xyj 

Similarly,  for  terms  of  degree  three,  we  obtain: 


r  /  \ 


#(2,2,3)  =  span  • 


V  0  / 


/  \ 

x2y 


V  0  / 


\x2y  J 


I  2\ 


V  0  / 

/  \ 


\xy  ) 


V  0  / 

/  \ 


V*3/ 


(   \  \ 


x3 


r  /  \ 


0 


5(2,2,3)  =  span  < 

A    °  / 
and  C{2, 2,3)  may  be  chosen  as 


x2y 


/  \ 

9 


v  0  / 


/  \ 


v0/ 


/  \ 


X3 

V x  J 


\xy ) 


\y  ) 


( (  \ 


C(2,2,3)  =  span  < 


IV  0  / 


/     \  1 


v  x2y  )  J 


For  k  >  4,  H(2,2,k),  i?(2,2,Jfc),  C(2,2,fc)  can  be  derived  accordingly.  Under  such  decom- 
positions, the  normal  form  up  to  degree  5  will  take  the  form: 


x    =    a2x2  +  a3x3  +  a4x4  +  a5x5 


y   =    -y  +  b2xy  +  b3x2y  +  b4x3y  +  b5x4y 
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The  actual  description  of  the  normal  form  is  not  unique.  Two  possible  descriptions  of  the 
normal  forms  up  to  degree  5  are  given  as  follows: 

x    =    -x3  -  4x5 

y   =    -y  +  4xy2  +  8xyi 

and 

x     -    -x3  -  Ax5 

y    -    -y  +  4xy2  +  20xy4. 

000 

The  fact  that  gk  €  C(n,  n,  k)  means  that  gk  is  just  a  linear  combination  of  the  basis 
elements  which  span  C(n,n,k).  To  determine  the  structure  for  gk,  we  need  to  choose  the 
space  C(n,n,k)  and  its  basis.  The  computation  shown  above  is  tedious,  prone  to  mistakes, 
and  not  systematic.  Fortunately,  one  can  easily  find  C(n,  n,  k)  and  the  associated  basis. 
Theorem  10  below  presents  a  method  for  finding  C(n,  n,  k)  for  a  given  matrix  A.  The  proof 
of  this  theorem  is  sketched  [4]. 

Theorem  10  [4]Ker{LkA.)  is  an  orthogonal  complementary  subspace  to  R(n,  n,k)  in  H(n,n,k), 
i.e., 

H(n,n,k)  =  R(n,n,k)®  Ker(LkA.), 

where  Ker(LkA.)  =  {{k  \  (LkA.£k)(y)  =  o}  and  A'  =  AT ,  i.e.,  transpose  of  conjugate  of  A. 

000 

Example  7  (Use  of  Ker(LA.) ) 
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We  use  Theorem  10  to  find  C(2, 2,  k)  for  systems  with  A 


Let  £*(»)  = 


/  \ 
0  1 


£2(2/1 .  2/2) 


.  Then 


i 


d(k2/dyi  dfi/dyt 
\ 

V\  dtf/dy2 
{  Vid&/dy»-&  J 


\ 


0  0 

V1  °J 


\  (  \ 


(  \ 
0  0 


V1  °/ 


£2(2/1,2/2) 


It  is  easy  to  see  that  the  solutions  to  the  equation  LkA.(Zk(y))  =  0,  i.e., 


2/i  dtf/dm  =  0 


are 


fi(yi,»a)    =  a*2/i 


^2(2/1^2)    =  *  2/2  +  Ml- 


r  / 


Thus  C(2, 2,  fc)  =  span  « 
form  up  to  degree  k  is 


I  V 


Vi 


vt1 2/2 


/     \  1 


Therefore  the  structure  of  the  normal 


Ui    =    2/2  +  ^22/?  +  ^32/?  +  •  •  •  +  ak  2/i 


2/2    =    (&22/?  + 022/1^2) +         + 022/^2) +  ---  +  (&*yi  +a*:2/i  12/2)- 
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ooo 

In  the  example  7,  it  is  observed  that  the  dimension  of  C(2,2,  k)  is  two  for  any  integer 
k  >  2.  However,  for  each  k  >  2,  there  are  three  monomials  of  degree  k  in  normal  forms. 
In  many  situations,  it  is  desirable  to  further  reduce  the  number  of  monomials  presented  in 
the  dynamics  of  normal  forms.  The  next  theorem  tells  how  one  can  make  such  reductions. 

Theorem  11  [4]Let  dk  =  dim  A'er(2^.)  and  {v^,- ■  ■  ,vdk]  be  a  basis  of  Ker(LkA,).  If 
{wi,  •  •  • ,  w<ik}  C  H(n,  n,  k)  satisfies 


Wj)  =  < 


Si  ±  0     ifi  =  j 
0  otherwise 


i,  j,  =  1,2, •••,<&, 


then  the  span{wu-  ■  ■  ,wdk)  C  H(n,n,k)  is  also  a  complementary  subspace  to  R(n,n,k)  in 
H{n,n,k).  000 


Example  8  (Further  Simplification  of  Normal  Forms) 


(  ( 


In  the  previous  example,  {vi,  v2}  =  < 
f  / 


\ 


Vi'1 2/2 


Take  {^1,^2}  =  < 


/       \  1 


V  *  1  J 


>  is  a  set  of  basis  of  C(2, 2,  A;). 


0 


yt1  V2 


0 


>.    It  is  easy  to  verify  that  for  i,  j  =  1,  2, 


(t;,-,  t0j)  ^  0  only  if  i  =  j,  using  the  definition  of  the  inner  product  introduced  in  [4]. 
Therefore,  a  different  normal  form  up  to  degree  k  is  given  by 


Vi    =  2/2 


(b2yl  +  022/12/2)  +  (b3yi  +  fl22/i22/2)  +  ---  +  (bkVi+aky1  1  yi)  ■ 


000 
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Example  9  (Convergence  of  Normal  Forms) 
Consider  the  system 

i\    =    — Xx  +  x\Xi 
x2    =  V2x2. 

One  of  the  full  normal  forms  is  given  by 

til    =  -2/1 
2/2    =  V2y2, 

which  is  obtained  using  the  transformation 

X2     =  2/2- 

Example  10  [10]  (Non-convergence  of  Normal  Forms) 

In  this  example  it  will  be  shown  that  the  normal  form  transformation  is  not  always 
convergent. 
The  system 

i\    =  x\ 

i2    =    x2  —  x\ 
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can  be  transformed  to 


Vi    =  V\ 


2/2     =  2/2 


using  the  transformation 


00 


X2     =     Hi  +  J>  "  ^  »*■ 


k=l 


Unfortunately,  this  transformation  is  not  convergent  in  any  small  neighborhood  of  the 


In  this  Section,  we  present  an  efficient  algorithm  for  deriving  normal  forms  and  the 
associated  transformations.  A  general  algorithm  for  obtaining  the  transformed  dynamical 
system  from  the  original  one  is  introduced  first. 
3.2.1    A  General  Algorithm  for  Coordinate  Transformation 

Consider  a  nonlinear  dynamical  system  depending  on  a  parameter  as  follows, 


origin. 


000 


3.2    Computation  of  Normal  Forms 


x  =  f(x,  a) 


(3.12) 


where  /  :  Rn  X  R  —  Rn  and  0  ^  a  €  R. 


Let 


x  =  u(y,a) 


(3.13) 
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be  a  transformation  (henceforth  called  "backward  transformation")  with  u(y,0)  =  y,  with 
its  inverse  (henceforth  called  "forward  transformation")  given  by 

y  =  v(x,a).  (3.14) 

Clearly,  we  also  have  v(x,0)  =  x. 

If  the  transformation  (3.13)  is  applied  to  the  dynamical  system  (3.12),  then  the  dynamics 
for  y  is 

.     v  A  dv(x,a)  .  (~  . 

y  =  g{y,a)  =     y       f(x,  a)|I=u(v,a).  (3.1o) 

It  is  desirable  to  have  an  efficient  algorithms  to  derive  equations  (3.15)  from  (3.12)  using 
(3.13)  and/or  and  (3.14).  To  simplify  the  description  of  the  algorithms  to  be  derived,  the 
following  partial  differential  equation  (PDE)  is  introduced: 

p-  =  U(u,a)Mv,*)  =  +  (3-16) 
oa 

Let  u(y,a)  be  the  (assumed  to  be  unique)  solution  to  the  PDE.  The  backward  transforma- 
tion is  then  constructed  as  x  =  u(y,a).  Later  it  will  be  shown  that  for  any  given  function 
u(y,a),  there  exists  a  unique  function  U(u,a)  satisfying  (3.16),  vice  versa.  The  algorithm 
for  computing  U(u,a)  from  a  given  u(y,a),  or  for  computing  u(y,a)  from  a  given  U(u,a) 
will  be  derived  later. 

In  practice,  only  the  first  few  terms  in  the  Taylor  series  of  /  (or  u,  v,  g,  U  respectively) 
in  a  are  important  in  stability  analysis.  However,  it  is  convenient  to  develop  the  method 
for  infinite  power  series  in  a  even  though  none  of  functions  /,  u,  v,  g,  and  U  will  necessarily 
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have  a  convergent  power  series  in  a.  This  will  present  no  difficulty  as  long  as  only  a  finite 
number  of  terms  of  the  series  are  calculated. 

Using  (3.16),  an  algorithm  for  deriving  (3.15)  from  (3.12)  is  given  in  the  following 
Theorem. 

Theorem  12  [15]  Let  f(x,a),  g{y,a)  and  U(u,a)  in  (3.12),  (3.15)  and  (3.16)  be  expanded 
as  follows: 


—  a 

f(x,a)   =    E/m(*)^  (3-17) 

m=0 

oo  m 

9(y,a)   =    Y.9m{y)a—V  (3-18) 

Um{u)a 

m=0 


m=0 

oo 


U(u,a)   =  (3-19) 


The  sequence        in  (3.18)  satisfies 


9i  =  /o'\  (3-20) 


where,  for  i  =  0, 1,2, ••  •,  the  sequence  fj™l  is  defined  recursively  via 


//0)    =    U  (3-21) 


i—m 


AH   =   /L^UEI  7    K/lr-!]'  m=l,2,...,i.  (3.22) 


j=0 


Here  (*-?")  =  jn^^l  "  irnomia/  coefficient  and  [/,  ff](i)  =  §*/  -  ff  5  is  tAe  lie 
bracket.  o  o  0 


This  algorithm  can  be  visualized  as  follows, 


fo  =  /f  =  90 

-> 

h  =  A(0) 

i 

1  _  AO) 

- 

A(1) 

-    fo  =92 

I 

i 

/3  -  J  3 

— ► 

r(2)         _      f(3)  _ 

1 

i  i 

-*     /a        ^  Vi 

1 

1  1 

The  dependency  between  {Um},  {fm},  {gm}  is  clear  from  this  graph. 
The  first  few  elements  in  the  sequence  {gm}  can  be  easily  written  out 

go    =  fo, 

gi    =    fi  +  [Uo,fo]  =  fi-[fo,Uo], 

52  =    f2  +  Wo,fi]  +  [Uo,  9i]-[foM, 

53  =    A  +  [Pb,  /a]  +  2  Pi ,  A]  +  [0o,  A  +  [^o,  A]  +  [0i ,  /o]  ] 

+[0i,5i]  +  [0o,52]-[/o,02]. 

In  general,  it  is  easy  to  verify  that 


<7m  =  rm  -  [/o,  Um-i] 
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where  Tm  only  depends  on  /0,  fi,  •  •  • ,  /m,  and  U0,  ■  •  • ,  Um-2-  The  similarity  between  (3.7) 
and  (3.23)  should  be  observed. 

The  next  Theorem  is  useful  for  deriving  an  explicit  expression  of  u(y,a)  when  the 
sequence  {Um}  is  known. 

Theorem  13  [15]  Assume  u(y,a)  is  the  solution  to  the  PDE  (3.16),  and 


—  G 

p(x,a)  =       pm(x)  —  . 


m=0 


(3.24) 


Let 


?(y,o)  =  p(u(y,a),a)  =  irniy)—;-) 


(3.25) 


m=0 


The  sequence  {qm}  in  (3.25)  satisfies 


(3.26) 


where,  for  i  =  0, 1, 2,  •  •  • ,  the  sequence  p\ 


is  defined  recursively  via: 


i—m 


(0) 
Pi  = 


Pi 


(3.27) 


n(m)      _       (m-1)  V- 


t—m 


j=0 


(3.28) 


Here  =  ^-ZZ-jY  is  the  binomial  coefficient  and  f  *  g{x)  =  ^f-g(x)  is  the  Lie 


derivative. 


ooo 
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Similarly,  the  algorithm  can  also  be  visualized  as  follows: 


(o) 

Po  =  Po   =  9o 


(o) 
Pi  =  Pi 


(o) 

P2  =  n 


(0) 
P3  =  P3 


i 

Po]  =  9i 

i 

P? 

f  1 

(2) 

— ►    Po    —  92 

i 

1 

(2) 

(3) 

-    Po    =  93 

1 

1 

1 

(2) 

(3) 

-  Pi 

i 

I 

i 

1 


The  dependency  between  {Um},  {pm},  {9m}  is  observed  clearly  from  this  graph. 
A  procedure  for  deriving  the  explicit  expression  of  u(y,  a)  is  given  as  follows. 

Procedure  1:  Applying  Theorem  13  with  p(x,a)  =  U(x,a)  yields  an  expansion  for  q(y,  a)  = 
U{u(y,a),a)  =  £m=o 9m(j/)^f-  Substituting  u(y,a)  and  U(u(y, a),a)  into  the  left 
and  right-hand  sides  of  the  PDE  (3.16),  respectively,  gives 

Wm+l  =  9m,  m  =  0, 1,2,3,  ••  •. 

000 

The  first  few  elements  in  the  sequence  {qm}  can  also  be  easily  written  out  as  follows: 


My)  =  io(y)  =  u0(y) 
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«*(»)   =   1l(y)  =  Ul(y)+^^Ut,(y) 


It  is  also  easy  to  verify  that  um  has  the  general  form: 

Um  =  tfm-i  +  Am,  m  =  1,2,3, (3.29) 

where  Am  only  depends  on  U0,  Ui,  •  ■  ■ ,  Um-2-  Therefore,  we  can  conclude  that  for  a  given 
sequence  {um},  there  exists  a  unique  sequence  {Um}  satisfying  the  PDE  (3.16),  vice-versa. 
Later,  we  will  need  the  following  Lemma  which  is  proved  in  [15]. 

Lemma  1  [15]  If  each  fmis  a  homogeneous  polynomial  map  of  degree  o/m  +  1,  and  each  Um 
is  a  homogeneous  polynomial  map  of  degree  ofm  +  2,  then  gm,  vm,  and  um  are  homogeneous 
polynomial  maps  of  degree  o/m+1.  000 

3.2.2    A  Procedure  for  Computation  of  Normal  Forms 

We  have  seen  that  the  normal  form  can  be  derived  using  a  sequence  of  near  identity 
transformations.  This  computation  method  is  very  inefficient  even  it  is  good  for  demon- 
stration of  the  concept  and  structure  of  normal  forms.  In  additional,  this  method  is  not 
appropriate  for  implementation  using  a  symbolic  computation  package  such  as  Maple  or 
Mathematica.  Alternatively,  the  algorithms  introduced  in  Theorem  12  and  13,  with  slight 
modifications,  can  be  used  to  derive  normal  forms  and  the  associated  transformations. 
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Assume  the  original  nonlinear  dynamical  system  is  given  by 


OO  1 

* = /w  =  £  /M+1(*)A,  **Rn-  (3-3°) 


m=0 


where  /^x)  =  Ax  with  A  =  I>/(0),  fm  £  #(n,  n,  m) .  Let  the  sequence  of  decompositions 


JT(n,  n,  m)  =  £(n,  n,  m)  ©  C(n,  n,  m) 


be  chosen.  We  need  to  obtain  the  transformations 


OO  -I 

x  =  u(y)=  5>m+1(y)-±j,  2/£i?n'  (3-31) 

m=0 


where  u1(y)  =  y  and  um  €  H(n,n,m),  as  well  as 


y  =  v(x)  =  £  ^m+1(x)^  (3-32) 

m=0 


where  ^(x)  =  x  and  vm  £  H(n,n,m),  so  that  the  dynamics  for  y  are  given  by: 


y  =  9(y)^d^f(X)Uu{y)-tr+\y)^  (3-33) 

m=0 


where  y*(y)  =  Ay,  and  gm  £  C(n,n,m).  First  of  all,  we  introduce  an  artificial  parameter 
by  scaling  the  variables  according  to  x  i — ►  ax,  y  \ — *  ay,  and  obtain 


x    =    /(ax)/ai/(x,a)=  £/m+1(x)— ,  (3-34) 

0 771. 

a  °°  am 

x    =    u(ay)/a  =  u(y,a)=  £  ttm+1(y)— ,  (3.35) 


m=0 
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y   =   v(ax)/a  =  v(x,a)=  £  vm+\x)  — ,  (3.36) 

0 771. 

y   =   g(ay)/a^g(y,a)=Y,Sm+\y)-^-  (3-37) 


m=0 


The  PDE  will  be,  according  Lemma  1, 


|H  =  U(u,a)  =  V;  tr+^t.)^,  «(y,0)  =  2/-  (3.38) 


For  the  purpose  of  computing  normal  forms,  v(-)  in  (3.32)  may  be  taken  to  be  the 
normal  form  transformation  and  the  dynamics  (3.33)  would  be  the  resulting  normal  form. 
Thus  the  process  of  computing  the  normal  form  dynamics  can  be  accomplished  recursively 
with  successive  computation  of  g'(-).  This  procedure  is  summarized  as  follows: 

Procedure  2:  Take  fm,  gm  and  Um-l  in  Theorem  12  to  be  fm+\  gm+\  and  Um+1  in  equa- 
tions (3.34),  (3.37),  and  (3.38),  respectively.  Beginning  at  Step  1,  with  m  =  1  and 
g1  =  f1,  each  pair  Um  and  gm  is  determined  iteratively.  At  Step  m,  a  parameter- 
ized form  for  Um  €  H(n,n,m)  is  assumed.  Note  that  g1,  ■  ■  ;  g"1'1  and  U2, 
Um~l  are  determined  from  previous  steps.  The  homogeneous  mapping  gm  is  then 
computed  using  the  algorithm  in  Theorem  12.  Restricting  gm  to  C(n,  n,  m)  will  fix 
some  coefficients  in  Um.  Other  coefficients  in  Um  are  free  to  be  chosen.  These  free 
coefficients  are  simply  set  to  zero  in  this  dissertation.  Once  all  of  the  coefficients  in 
Um  are  determined,  both  um  and  gm  are  also  determined.  OOO 

3.3    Further  Discussion  about  Normal  Forms 

We  have  seen  that  (with  f0(x)  =  Ax,  gm  =  gm+1,  Um-\  =  Um+1) 

gm  =  rm  +  [um-i,  fo]  =  Tm  -  [/o,  um-i] 
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=    rm-I?+1tfm_i>  (3-39) 

where  Tm  is  the  terms  of  degree  m  +  1  after  transformation  x  =  YlT=o  Uk(y)/M  has  been 
applied  to  the  original  nonlinear  dynamical  system.  For  the  linear  operator  the  space 

H(n,n,m+  1)  has  the  following  two  decompositions  (one  in  the  domain,  the  other  in  the 
range): 

ff(n,n,ro+l)    =    H?(n,n,m+1)Q  H?(n,n,m+ 1) 
=    Hp(n,n,m  +  l)§  H?(n,n,m+  1), 

where  H? (n,  n,  m  + 1)  is  the  null  space  of        ,  fff  (n,  n,  m  + 1)  is  the  range  of         ,  i.e., 

JTftn.i^m+l)  =  R{L^)  =  R{n,  n,  m  +  1), 

fff(n,n,m  +  l)  =  iV(X^.+1)  =  C*(n,  n,  m  +  1), 

#f(n,n,m+l)  = 

JT?(n,n,m+l)  =  iV(X^+1). 

Let  Pj0  be  the  projection  of  .ff(n,  n,  m  +  1)  onto  /rf(n,  n,  m  +  1)  along  #f  (n,  n,  m  + 
1).  Similarly,  let  Pf1  be  the  projections  of  H(n,n,m+  1)  onto  H^(n,n,m  +  1)  along 
^2fi(n,n,m+l).  Let  be  the  inverse  of  L%+1  taking  #f(n,  n,  m+1)  to  #P(n,  n,  m+ 

1).  A  rigorous  definition  of  the  normal  form  is  given  as  follows. 

Definition  i  [15]  Let  (3.30)  be  the  original  nonlinear  dynamical  system.  Equation  (3.33)  is 
a  normal  form  for  (3.30)  relative  to  the  matrix  A  and  projections  Pf* ,  PXH  if  the  elements 
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gm,  Um  satisfy 

gm(y)    =  (/-Jf)rro(y) 
Um-i(u)    =  (^+1)"1^r"l(i') 

where  Tm  is  defined  in  the  equation  (3.39).  000 

Since  :  H(n,n,m)  -»  H(n,n,m)  is  a  linear  operator,  it  can  be  represented  using  a 
matrix  of  appropriate  dimension.  It  is  also  desirable  to  know  the  eigenvalues  of  this 
operator.  Remember  that  a  complex  number  A  is  a  eigenvalue  of  L%  if  and  only  if  there 
exists  a  nonzero  Pm  €  H{n,n,m)  (possibly  over  complex)  such  that  L%Pm(x)  =  APm(x). 

The  following  theorems  are  proved  in  [15]. 

Theorem  li  [15]  Let  an  n  x  n  matrix  T  be  invertible  and  B  =  T~XAT.  The  eigenvalues  of 
the  linear  operators  U\  and  Lg  are  identical,  i.e., 

In  other  words,  the  spectrum  of  the  linear  operator  L™  is  invariant  under  a  linear  change 
of  coordinates. 

Theorem  15  [15]  Let  the  eigenvalues  of  a  given  nx  n  matrix  A  be  {Ai,  A2,  •  •  • ,  An},  then 


*(La)  =  IE  A<°'  -  A,  |j  =  1,  2, 3,  •  •  -  n;  f>t  =  m;  a,  >  0  i 
li=i  «=i  > 


(3.40) 


000 
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We  now  study  how  to  derive  the  matrix  Z^.  H(n,  n,  m)  is  a  linear  space.  The  dimension 

of  H(n,  n,  m)  is  equal  to  n^&frffi  ■  The  basis  for  ^  (n' n' m)  can  be  determined  as  follows. 
As  we  have  seen  before,  the  set  of  monomials  {xa  |  |a|  =  m}  is  a  basis  for  H(n,  1,  m),  and 
dimension  of  H(n,  l,m)  is  (n+™_1).  Let  {eu  ■  •-,«„}  be  the  natural  basis  for  Rn,  i.e.  the 
j"1  entry  of  the  column  vector  ej  is  one,  while  all  other  entries  are  zero.  A  basis  for  the 
space  H(n,  n,  m)  is  given  by 

{xaeJ\\a\  =  m,l<j<n},  (3.41) 

with  dm  =  dim  H(n,  n,m)  =  n-  (n+^_1).  We  order  the  basis  elements  (described  in  (3.41)) 
of  H(n,  n,  m)  in  the  reverse  lexicographic  order,  i.e.,  iae,  proceeds  x3€j  if  and  only  if  the 
first  nonzero  entry  in  the  row  vector  [i-j,  cti  -  ft,  •  •  • ,  an  -  /?„]  is  positive.  Therefore, 
a  correspondence  between  the  basis  elements  of  Rdm  and  H(n,n,m)  can  be  established. 
When  both  n  and  m  are  big,  finding  and  ordering  the  basis  for  H(n,n,m)  is  not  a  easy 
job.  Fortunately,  we  can  use  symbolic  computation  software  such  as  Maple  to  generate  the 
basis. 

Example  11  (Derivation  of  Basis) 

The  basis  for  H(2, 1, 2)  is  {x2,  xy,  y2}  in  reverse  lexicographic  order  and  dim  #(2, 1,2)  =  3. 
The  basis  for  H '(2,2,2)  in  the  reverse  lexicographic  order  is 


and  dim  H(2, 2, 2)  =  6. 


OOO 
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Once  the  basis  for  a  linear  space  is  fixed,  the  matrix  representation  of  the  linear  operator 


I?  can  be  derived  by  applying  it  to  each  basis  element  of  H(n,  n,  m). 


Example  12  (Derivation  of  the  Matrix  Representation  of  the  Linear  Operators  L%) 


For  A  = 


0  -1 


1  0 


with  a  pair  of  critical  eigenvalues,  ±j,  the  matrix  representation  for 


/ 


the  operators  L2A,  LA  can  be  expressed  by 


/  \ 
0  10-100 


-  2  0    2    0     -1  0 

0  -1    0     0      0  -1 

1  0  0  0  1  0 
0  10-202 
0  0  10-10 


,  and 


V 


0 

1 

0 

0 

-1 

0 

0 

0 

-3 

0 

2 

0 

0 

-1 

0 

0 

0 

-2 

0 

3 

0 

0 

-1 

0 

0 

0 

-1 

0 

0 

0 

0 

-1 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

-3 

0 

2 

0 

0 

0 

1 

0 

0 

-2 

0 

3 

0 

0 

0 

1 

0 

0 

-1 

0 

/ 


The  resulting  eigenvalue  sets  are  {3  j,  -3  j,  j, -j}  and  {0, 0,4  j,  -4  j,  2  j,  -2  j, 2  j,  -2j} 
respectively. 

3.4  Summary 

In  this  chapter,  normal  form  theory  was  introduced.  The  concept  and  structure  of 
normal  forms  were  first  illustrated  using  near  identity  change  of  coordinates.  An  efficient 
algorithm  for  the  computation  of  the  normal  form  was  presented.  The  problem  of  simplifi- 
cation of  nonlinear  dynamics  could  be  reduced  to  the  study  of  the  linear  operator  LA. 


CHAPTER  4 

RELATIONSHIP  BETWEEN  CENTER  MANIFOLD  AND  NORMAL  FORM  THEORIES 

As  we  will  see,  center  manifold  and  normal  form  methods  are  two  major  model  reduction 
techniques  in  study  of  critical  nonlinear  systems.  In  this  chapter,  the  relationship  between 
these  two  methods  will  be  explored.  As  we  already  knew,  center  manifolds  for  a  given 
critical  nonlinear  system  are  not  unique.  However,  it  will  be  proved  in  Section  4.1  that 
all  center  manifolds  for  a  given  critical  nonlinear  system  are  equivalent  in  some  sense. 
In  Section  4.2,  normal  forms  of  critical  nonlinear  systems  are  studied.  The  relationship 
between  center  manifold  and  normal  form  theories  is  presented  in  Section  4.3.  Section  4.4 
includes  an  example  to  demonstrate  the  equivalence  of  center  manifold,  normal  forms  of 
critical  nonlinear  systems,  and  the  relationship  between  center  manifold  and  normal  form 
theories.  A  brief  summary  in  section  4.5  concludes  this  chapter. 

4.1    Equivalence  of  Center  Manifolds 

As  we  mentioned  before,  center  manifolds  need  not  to  be  unique  for  a  given  nonlinear 
system.  In  this  section,  we  will  prove  that  all  center  manifolds  for  a  given  Ck  nonlinear 
system  are  A; -jet  equivalent.  First  of  all,  the  concept  of  k-jet  equivalence  is  defined. 

Definition  5  [17]  Let  f,  g  :  D  C  Rn  —>  Rm  be  two  Ck  maps  with  k  >  1,  being  an  integer. 
The  maps  f  and  g  are  k-jet  equivalent  on  a  setQ.  C  D  if  and  only  if  their  partial  derivatives 
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up  to  order  k  coincide  on  D: 


for  all  x0  £  D,  and  non-negative  integers  [h,  k2,  ■  ■  • ,  kn],  with  0  <  E"=i  fc<  -  k-  In  the 
other  words,  the  Taylor  series  expansions  of  f  and  g  on  D  coincide  up  to  degree  k.  000 

The  following  lemma  is  needed  later.  The  proof  of  this  lemma  is  found  in  [15]. 

Lemma  2  [15]  For  a  given  nxn  matrix  A  and  a  given  m  x  m  matrix  B,  with  eigenvalues 
{Ai,  •••,*„},  and{m,---,nm},  respectively,  define  the  linear  operator  UAB  :  H{m,n,r)  - 
H(m,n,r)  via 

L\B{Pr{*))=d-^1Bx-AF{x),  (4.1) 
where  Pr(x)  £  H(m,n,r).  Then  the  spectrum  of  UA  B  is  given  by 


°(La,b)  =  {  E  »W  -  Aj  |  j  =  1,2,  •  •  -  n;  J>  =  '  with  qi>o\.  (4.2) 
I  i=i  <=i  > 


000 

Consider  a  critical  Ck  nonlinear  system  defined  on  some  neighborhood  Q.  of  the  origin 
of  Rn,  (k  >  2).  Without  loss  of  generality,  this  system  may  be  described  as  follows; 

xc    —  Acxc fc(xc,  xs) 

=  Acxc  +  f2c{xc,  *,)  +  ..-+  /cfc(xc,  x,)  +  0(|*|*+1),  (4.3) 

%s    ~  A  x3  -f-  fs{xc,x3) 

=  A*xs  +  fi(xe,  *,)  +  •■•  +  /*(xc,  x.)  +  0(|x|*+1),  (4.4) 
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where  xc  €  Rnc,  z,  e  Rn',  nc  +  ns  =  n,  /'  €  H(n,  nc,  t),  and  /j  €  #(n,  n„  »')■  From  center 
manifold  theory,  there  exists  a  Ck  center  manifold, 

xs  =  h(xc)  =  h\xc)  +  /i3(xc)  +  •  •  •  +  fcfc(*e)  +  0(Nfc+1),  (4.5) 

for  the  system  (4.3),  (4.4)  on  some  neighborhood  ft'  C  Anc,  such  that  the  local  stability  of 
the  system  (4.3),  (4.4)  is  equivalent  to  the  local  stability  of  the  subsystem 

xe=  Acxc  +  fc(xc,h(xc)).  (4.6) 

Theorem  16  All  center  manifolds  for  (4.3),  (4.4)  are  k-jet  equivalent  on      .  OOO 

Proof:  A  parameterized  function  xs  =  h(xc)  =  h2(xc)  +  h3{xc)  +  •  •  ■  +  hk(xc)  +  0(|*c|*+1) 
is  assumed  first.  If  xs  =  h(xc)  is  a  center  manifold  for  (4.3),  (4.4),  by  definition,  it  satisfies 
following  Partial  Differential  Equation  (PDE), 


(Acxc  +  fc(xc,  h(xe)))  =  Ash(xc)  +  /.(*„  h(xc)).  (4.7) 

oxc 


More  specifically, 


d{h\xc)  +  h*{Xc)  +  ■■■)  (AcXc  +  h(Xe))  +  /c3(Xc)  h(Xc))  +\  (48) 

oxc  v  ' 

=  A°h(xc)  +  fl(xc,  h(xc))  +  /s3(xc,  h(xc))  +  ■■-.  (4.9) 


It  is  clear  that  the  homogeneous  terms  of  degree  r  <  k  on  the  left-hand  side  of  the 
equation  (4.9)  only  depend  on  h2,  •  •  •,  hT  and  Ac,  /c2,  •  •  • ,  //.  Similarly,  the  homogeneous 
terms  of  degree  r  <  k  on  the  right-hand  side  of  the  equation  (4.9)  only  depends  on  h2,  •  •  • ,  hr 
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and  Aa,  /2,  •  •  • ,  //.  With  a  little  work,  we  can  show  that  hr  satisfies  the  following  PDE. 


m^-Acxc  =  A°hr(xc)  +  IT(*C)  (4.10) 

OX  c 


for  some  function  IT(xc)  which  depends  on  h2,  ■  ■  • ,  hr~l,  /c2,---,/cr  and  /2,  •••,/;.  Notice 
that  all  eigenvalues  of  Ac  have  zero  real  parts,  while  all  of  the  eigenvalues  of  As  have 
negative  real  parts.  From  the  Lemma  2,  the  linear  operator  UA,  Ac  does  not  have  any  zero 
eigenvalues.  This  means  that 


LTA.,Ahr{xc))  =  ^jf^Acxc  -  A'hr(xc)  =  r (xe),  (4.11) 


has  an  unique  solution  hr(xc)  for  any  function  IIr(xc)  €  H{nc,ns,  r)  where  r  =  2, 3, 

<>oo 

Example  13  (Equivalence  of  Center  Manifolds) 
We  have  shown  that  for  any  real  number  cl5C2, 


x„  =  h(xc,ci,c2)  =  < 


Ciexp(-5XC2),   xc  >  0 
0  xc  =  0 

c2exp(-lx72),   xc  <  0 


is  a  center  manifold  for  the  system 


xc  =  -xl  (4.12) 
i,    =    -x,.  (4.13) 
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It  is  very  straightforward  to  show  that  hk  (xc,c1,c2)  =  0  for  any  integer  k  >  0. 

4.2    Normal  Forms  of  Critical  Nonlinear  Systems 

Among  other  applications,  normal  forms  are  capable  of  bringing  out  the  essential  fea- 
tures of  nonlinear  systems.  In  this  section,  we  use  normal  form  technique  to  simplify 
description  of  critical  nonlinear  systems.  As  we  proceed,  the  reason  the  normal  form  tech- 
nique is  called  a  model  reduction  method  will  become  evident. 

First,  let  us  consider  an  analytic  critical  nonlinear  systems  over  complex  of  form 

x=Ax  +  f(x)  =  Ax  +  f2(x)  +  f{x)  +  ---,    xeCT.  (4.14) 

Definition  6  [4]  The  eigenvalues  {Al5  A2,  •  •  • ,  An}  of  A  are  resonant  if,  for  at  least  one  value 
of  r,  1  <  r  <  n, 

n 

holds  for  some  non-negative  integers  qd=       •  •  • ,  qn]  satisfying  Ya=i  Qi  >  2- 

Remark  6  The  eigenvalues  of  any  real  matrix  with  critical  eigenvalues  are  necessarily  res- 
onant. 

Let  e  =  {ei,  •••en)  be  the  standard  basis  for  Rn,  so  that  the  ith  entry  of  et-  is  1  while 
all  other  entries  are  zero.  Let  xq  =  arf  x\?  •  ■  •  xqnn  be  a  scalar  monomial. 

Definition  7  [4]  The  vector  monomial  xqej  is  resonant  with  respect  to  an  upper  triangular 
matrix  A  (possibly  complex)  with  (ordered)  eigenvalues  X\, . . .,  An,  if,  for  those  non-negative 
integers  [q\, ... ,  qn]  indexing  xq , 

Aj-f>A,  =  0.  (4.15) 
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The  integer  r  =  Yli=i  9«>  >  2  is  called  the  degree  of  the  resonant  monomial.  OOO 

The  following  two  lemmas  are  proved  in  [4]. 

Lemma  3  [4]  If  A  is  an  upper  triangular  matrix  with  diagonal  elements  {Ai,  A2,  •  •  • ,  A„}, 
then  a  normal  form  up  to  degree  k  for  the  system  (4-H)  is  given  by 

j=  Ax  +  f{x)  +  ■■■  +  f(x)  +  0(|x|r+1)  (4.16) 

where  TO*))  *  =  2, is  made  up  of  linear  combinations  of  all  resonant  monomials. 

OOO 

Remark  7  The  actual  description  ofT(x)  is  not  unique.  It  depends  on  A  and  other  factors. 
On  the  other  hand,  not  all  resonant  monomials  of  degree  i  are  included  in  f'(x).  OOO 

Lemma  I  [4]  Assuming  that  the  following  equation 

3F=  Ax  +  f{x)  +  ■■■  +  f(x)  +  0  (\x\k+1)  (4.17) 

represent  a  normal  form  up  to  degree  k  with  respect  to  A.  Then  after  a  linear  change  of 
coordinates  x  =  Ty,  the  resulting  equation 

y=  By  +  f(y)  +  ■■■  +  g\y)  +  0  (\y\k+1)  ,  (4.18) 

will  be  a  normal  form  up  to  degree  k  with  respect  to  B  :=  T_1  AT ',  where  f{y)  =  T-1/*  (?!/), 
£  =  2,  -  -  - ,  *r.  OOO 

Remark  8  In  this  dissertation,  only  real  systems  are  considered.  Some  real  matrices  have 
complex  eigenvalues  so  that  their  Jordan  forms  are  complex  matrices.  After  deriving  the 
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structure  of  normal  forms  over  the  field  of  complex  number,  the  structure  over  the  field  of 
real  numbers  will  be  can  be  determined  easily.  OOO 

Consider  a  real  critical  nonlinear  system  as  follows: 


xr  = 


x,  = 


Acxc  +  fc(xc,xs) 

Acxc  +  f2c{xc,  xs)  +  /c3(xc,  x.)  +  •  •  • ,   xc  G  Rnc 
Asxs  +  f3{xc,xa) 

Asxs  +  f2a{xc,  x.)  +  f*(xc,  *.)  +  •••,    xs  e  Rn' . 


(4.19) 


(4.20) 


The  following  theorem  describes  the  structure  of  normal  forms  for  the  system  (4.19)-(4.20). 

Theorem  17  All  normal  forms  up  to  degree  k  for  the  system  (4.19)-(4.20)  have  the  following 
structure; 


xr  ■ 


x,  = 


(4.21) 


Acxc  +  fc(xc,xs) 

A°Xc  +  Pc{Xc)  +  fl(Xc)  +■■■  +  fc(*c)  +  0(\XC,  Xs\k+l) 

A"xs  +  £(xc, x,)  +  j*(ie,  *,)+>  ••+  j*(Se)  +  0(\xc,  xs\k+1 )  (4.22) 


with  fl(xc,xs  =  0)  for  i  =  2,3,4, 

/ 

Proof  of  Theorem  17:  Let  A  = 


matrix  T  = 


Tn  0 

o  r22 


\ 


Ac  0 
0  A" 


\ 

.  There  exists  an  invertible  (possibly  complex) 


such  that  A  = 


/  \ 

/ 

Ac  0 

=  T~XAT  = 

^  0  A'; 

iffiWu  0 
0  igA'Tn 
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where  Ac  and  As  are  upper  triangular  matrices  with  (possibly  complex)  diagonal  elements 
being  {AC)1,  •  •  • ,  AC)7lc}  and  {As>1,  •  •  • ,  As,nj},  respectively.  After  the  linear  transformation 
x  =  Ty  is  applied  to  the  system  (4.19),  (4.20),  the  dynamics  for  y  become 

yc  =  Acyc  +  gc(yc,ys) 

=  Acyc  +  g2c(y^ys)  +  9c(yc,ys)  +  ---,  (4.23) 

ys  =  Asys+gs(yc,ys) 

=  Asys+g2s(yc,ys)  +  g3s(yc,y  *)  +  ■■••  (4-24) 

Assume  that  a  normal  form  for  the  system  (4.23)-(4.24)  (up  to  order  k)  has  the  form: 

yc  =  Acyc  +  gc(yc,ys) 

=  Acyc  +  fc(yc,ys)  +  fc(yc,ys)  +  ---  (4-25) 

Vs  =  ^3y3  +  9s(y^ys) 

=   Asys  +  fs(yc,ys)  +  fs(yc,ys)  +  ---.  (4-26) 
From  Lemma  3,  gkc(yc,  2/s)  only  contains  monomials 

n 

y    —  2/c,l         yc,nc  f s,l  f»,n,  >  Z-f  ** 

1=1 


such  that 


Kj  =  E  «A<*  +  £  9nc+, A,,t- ,  1  <  j  <  nc.  (4.27) 

i=l  t=l 

Recall  that  the  real  parts  of  Ac,,  (or  A,,,-)  are  zero  (or  negatively  respectively).  To  make 
equation  (4.27)  hold,  it  is  necessary  to  have  qUc+i  =  •  •  •  =  qn  =  0,  hence  gk(yc,  y,)  =  <7c(!/c)- 
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Similarly,  g^Vclls)  only  contains  monomials 

n 

,.q  _  „9J    .  .  .  „9nc  „9"c+l  .  .  .  yqn  „.  _  U 

V    -  Vc.l         Vc,nc  Vs,l  Vs,ns  I 

1=1 

such  that 

ASJ  =  ^ftAe.i  +  X>„c+,ASii ,  1  <  j  <  nt.  (4.28) 
i=i  i=i 

Therefore,  it  is  necessary  for  E"=nc+i  9«'  ^  0, hence  gks(yc,ys  =  0)  =  0.  The  linear  trans- 
formation x  =  Ty  applied  to  (4.25),  (4.26)  will  yield  the  system  (4.21),  (4.22).  OOO 

4.3  Relationship 

The  main  results  of  this  section  show  the  relationship  between  center  manifold  and 
normal  form  theory.  We  first  derive  center  manifolds  for  the  system  (4.21),  (4.22). 

Theorem  18  Any  center  manifold  for  the  system  (4-21),  (4.22)  satisfies  xs  =  0(|xc|fe+1). 

   2  — 3   

In  other  words,  the  k-jet  equivalent  of  any  center  manifold  h(xc)  =  h  (xc)  +  h  (xc)  H  h 

t(xc)  +  0(|xc|*+1)  must  satisfy  V{xe)  =  0,  for  i  =  2,  . . . ,  k.  OOO 

The  proof  of  Theorem  18  is  very  similar  to  that  of  Theorem  16,  except  in  this  case 
n'(-)  =  0,  for  i  =  2,  in  the  PDE  (4.7).  With  this  theorem,  the  dynamics  for  the 
system  (4.21),  (4.22)  along  any  center  manifold  x,  =  h(xc),  up  to  order  k,  satisfy 

jc    =    Acxc  +  Jc(xc,h(xc)) 
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Therefore,  the  stability  of  the  system  (4.21),  (4.22)  is  determined  by  the  stability  of  the 
subsystem  (4.21).  That  is  the  exact  reason  the  normal  form  method  is  called  a  model 
reduction  technique. 

The  relationship  between  the  two  subsystems  (4.6)  and  (4.21)  (up  to  degree  k)  is  given 
in  the  following  Theorem. 

Theorem  19  One  of  the  approximate  normal  forms  of  the  subsystem  (4.6),  (up  to  order  k), 
is  given  by  (4.21).  000 

To  derive  a  normal  form  up  to  degree  fc,  the  transformation  of  form 

xc   =    xc  +  u2c(xc,xs)  +  ---+ukc{xc,xs)  (4.29) 


X  •    —  x_ 


+  u2s(xc,xs)  +  ---  +  uk(xc,xs).  (4.30) 


is  sufficient.  The  center  manifolds  for  the  system  ,  as  well  as  the  transformation  between 
(4.6)  and  (4.21)  can  be  derived  from  equations  (4.29),  (4.30). 

Theorem  20  The  transformation  xc  =  xc  +  i>2(xc)  +  ■■■  +  i>k{xc)  taking  (4.6)  to  (4-21)  is 
obtained  by  setting  xs  =  0  in  equation  (4-29).  Any  center  manifold  xs  =  h(xc)  (up  to  degree 
k)  can  be  derived  from  equations  (4-29),  (4-30),  by  setting  xa  =  0.  OOO 


4.4    A  Example 

Example  U  (A  CASE  2  system) 
Consider  following  system: 

x    =    -y  +  xz2  +  xyz  +  x2  +  z2,  (4.31) 

y    =     I  +  22  +  I2  +  ;rz25  (4>32) 
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z   =    -z  +  xy  +  z  , 


(4.33) 


with  A  = 


(4.33)  is 


f  0    -1  0 


1     0  0 


vo    0  -1, 


.  The  center  manifold  (up  to  degree  3)  for  the  system  (4.31)- 


z  =  h(x,y)  =  -^  +  if+V4  +  T  +  °(W*>V\\4)- 
o       5       5  5 


The  dynamics  (up  to  degree  3)  along  the  center  manifold  is 


(4.34) 


x  =  -2/  +  :r2  +  0(||z,2/||4) 
y   =   x  +  x2  +  0(\\x,y\\4). 


(4.35) 
(4.36) 


The  normal  form  up  to  degree  3  for  the  system  (4.35)-(4.36)  is  given  by 


x  = 


y  - 


-y-l(x2  +  2/2)  (31 -7f) 
x-±-(x2  +  f)(7x  +  3y) 


z    -  -z. 


(4.37) 
(4.38) 
(4.39) 


The  center  manifold  for  the  system  (4.37)-(4.39)  is,  up  to  degree  3, 


2  =  0 


(4.40) 


The  backward  transformations  which  transforms  (4.31)-(4.33)  to  (4.37)-(4.39)  is 


y  = 


z  = 


6  y2z     4  z3 

+  — —  +  

25  5 


yz 

20 


x2y     43  xyz 


2f 
3 


— 2 

xz* 


6  x2z 
'  To"      25  + 


x2 


2xy  2f 


21  xyz  x3  xz' 
~ 50       T  +  ~20 


4 

2xy2 
9 

z22z3 
5  5 
Ufz 


50 


7,3 


13  y 
36 


5x2y     5xy2  + 


50 
llx2z 


18 


12 


_   *2    *y  ,  #2    -2  , 


U22     2  xz2     2 1/2^ 

+  -  h  +  — — 

T  5        5  5 

2f_  _^y  4xf 

+  15       15        5  5 


50 
2  x3 

2  xi/2 


x2z 


The  forward  transformation  is 


x2     2sy    2g2     3z2  23s3 
I_X+3        3         3        5  36 

5x2y     12x2z     2xy2     31  xyz  xz 
~~L2~  +    25    +_9~+    50    +  6 
5x2y     12x2z    2sy2     31  xyz  xz 
12~  +    25         9         50  6 
x2     2xy     2  y2     z2  x3 

y  =  y-j-—-~r+j+T 

25x2y     lx2z     11  xy2      17 xyz 

+  _ 18         50~+     12     +  50 

7  xz2     29  y3     7y2z  2yz2 

12   +   36        50  3 

i2     xy     y2      2     x3  ^ 

z    =    z+   +  z  -  —  +  — 

5       o       5  5  o 

xz2     4y2z     2yz2  ^ 
 z    —     -  -     "T  z 
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The  backward  transformation  which  brings  (4.35)-(4.36)  to  (4.37)-(4.38)  is  obtained 
from  equations  (4.41)-(4.42)  by  simply  setting  2  =  0,  that  is 

7x3     x2y     2y3  2xf 
~36~~   4   ~   3   +  9 
x3     13  y3     5x2y  5xy2 
T  ~    36   ~    18    +  12 

The  forward  transformation  which  brings  (4.35)-(4.36)  to  (4.37)-(4.38)  can  be  derived 
in  many  different  ways.  First,  the  algorithm  in  the  Chapter  3  can  be  used.  Alternatively, 
it  can  be  derived  by  solving  equations  (4.44)-(4.46)  together  with  z  =  0.  A  small  trick  can 
be  used  as  follows.  From  equation  (4.46),  we  have 

,mJ*a.U,*+***j£-&!.+A.  (4.47, 

y  5       5       5  555         5  5  J 

Next,  equation  (4.47)  is  substituted  recursively  into  equations  (4.44)-(4.45)  until  z  is  only 
occurs  in  terms  of  degree  4  or  higher.  After  truncating  terms  of  degree  4  or  higher,  we 
obtain  the  forward  transformation  as  follows: 

x2     2xy     2y2  x2y     5y^     x^_     13  xy2 

T~T"  +  ~3  6        IT  +  18  ~  36 

2xy     2y2     x2  29x2y     2y3     2xy2  x3 
3        3       3        36         9         3  12' 


x  = 


y  = 


x'  2xy 


3 

2xy 


x 

V+  —  + 
"     3  3 


+ 


2f_ 
3 

2? 
3 


The  center  manifold  (4.34)  (up  to  degree  3)  can  be  solved  from  equations  (4.41)-(4.43) 
by  simply  setting  z  =  0,  or  derived  through  recursive  substitution  using  the  equation  (4.47). 


4.5  Summary 

In  this  chapter,  the  equivalence  of  center  manifolds  for  a  given  critical  nonlinear  system  is 
studied  first.  Next,  the  structure  of  normal  forms  for  critical  nonlinear  systems  is  derived. 
The  reason  the  normal  form  method  is  called  a  model  reduction  technique  is  explained. 
The  relationship  between  center  manifold  and  normal  form  theories  is  explored.  Finally,  a 
example  is  provided  to  illustrate  this  relationship. 


CHAPTER  5 

CONSTRUCTION  OF  LYAPUNOV  FUNCTIONS  -  I 

This  chapter  is  concerned  with  the  construction  of  Lyapunov  functions  of  critical  non- 
linear dynamical  systems  if  they  are  LAS.  As  we  mentioned  before,  a  necessary  condition 
for  a  given  nonlinear  dynamical  system  x  =  f(x)  to  be  LAS  is  that  the  eigenvalues  of  its  Ja- 
cobian  matrix  A  =  Df(0)  are  either  stable  or  critical.  If  the  Jacobian  matrix  possesses  only 
stable  eigenvalues,  then  a  Lyapunov  function  for  the  corresponding  linear  system  x  =  Ax 
is  also  a  Lyapunov  function  for  the  nonlinear  system  itself.  The  only  situation  which  needs 
further  investigation  is  when  the  Jacobian  matrix  possesses  stable  eigenvalues  and  at  least 
one  critical  eigenvalue.  Therefore,  only  critical  nonlinear  systems  will  be  considered  in  this 
chapter.  In  Section  5.1  we  will  present  a  procedure  for  constructing  Lyapunov  functions 
for  general  critical  nonlinear  systems.  To  apply  this  procedure,  it  is  required  that  there 
exist  a  full  normal  form  for  the  given  critical  nonlinear  system.  In  Section  5.2  an  explicit 
algorithm  for  construction  of  Lyapunov  functions  for  nonlinear  systems  in  two  special  cases 
is  provided.  This  algorithm  only  requires  an  approximate  normal  form.  Examples  are  given 
in  Section  5.4,  to  demonstrate  use  of  the  algorithm  along  with  some  discussions.  The  major 
points  in  this  chapter  are  summarized  in  Section  5.5. 

5.1    General  Procedure  For  Constructing  Lvapunov  Functions 
Consider  a  critical  nonlinear  system  described  below: 

xc   =    A°xc  +  fc(xc,xa) 
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=    Acxc  +  f?(xc,xs)  +  f?(xc,xs)  +  ---,  xc<ERn<,  (5.1) 
x3    —    A  xs  -\-  /s(xc,ig) 

=    Asxs  +  tf(xc,xs)  +  f*(xc,xs)  +  ---,   xseRn>.  (5.2) 

If  a  full  normal  form  for  the  system  (5.1)  and  (5.2)  exists,  the  normal  form  will  have  the 
following  structure: 

xe    =    Acxc  +  fc(xc) 

=  Aexc  +  f?(zc)  +  f!(xe)  +  '-;  (5.3) 
2,    =    Asxs  +  f,(xe,xa) 

=    Asxs  +  f?(xc,xs)  +  f?(xc,xs)  +  ---,  (5-4) 

where  for  i  =  2, 3, 4,  •  •  •,  /j'(xc,  xs  =  0)  =  0. 

In  this  section,  a  general  method  for  the  construction  of  Lyapunov  functions  for  systems 
of  the  form  (5.3),  (5.4)  which  are  known  to  be  LAS  will  be  provided.  This  construction 
depends  on  Lyapunov  functions  Vc(xc)  and  Vs(xs)  associated  with  the  critical  and  stable 
dynamics,  respectively.  The  following  lemma  shows  the  stable  dynamics  (5.4)  are  always 
LAS  for  ||xc||  small  enough. 

Lemma  5  If\\xc\\  is  small  enough,  the  dynamics  of  the  stable  states  (5.4)  are  LAS.  OOO 

Proof:  Since  7c(xc,  0)  =  0  for  i  =  2, 3, 4,  •  •  •  in  (5.4),  the  dynamics  for  xs  can  be  rewritten 
as 

i4=  Asxs  +  G(xc,xs)xs  (5.5) 
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where  G(x~c,xs)  is  an  ns  x  ns  matrix  whose  entries  are  continuous  functions  of  the  com- 
ponents of  xc  and  xs,  with  G(0,0)  =  0.  Let  P  =  PT  >  0  be  a  real  matrix  such  that 
(AS)TP  +  PA'  <  0.  Taking  V„(xs)  =  x^Pxs  as  a  Lyapunov  function  candidate  for  (5.4) 
yields 

Ws(xc,xs)  =}Vs(xs)   =   xT{(As)TP  +  PAs}xs  +  xJ{GT(xc,xs)P 

+PG(xc,xs)}xs.  (5.6) 

There  exists  some  neighborhood  of  xs  =  0,  such  that  in  this  neighborhood,  each  entry  of  the 
matrix  GT(xc,xs)P  +  PG(xc,xs)  will  be  arbitrarily  small,  if  ||xc||  is  small  enough.  Thus, 
GT{xc,xs)P  +  PG(xc,xa)  can  be  treated  as  a  small  perturbation  to  the  matrix  (AS)T P  + 
PA",  which  remains  negative  definite  under  sufficiently  small  perturbations.  Therefore,  for 
||xc||  small  enough,  Ws(xc,0)  =  0  and  Ws{xc,xs)  <  0  for  xs  /  0  in  some  neighborhood  of 
xs  =  0.  Thus,  the  subsystem  (5.4)  is  LAS  for  ||xc||  small  enough.  Furthermore,  V,(xs)  = 
xJPxs  is  a  Lyapunov  function  of  the  subsystem  (5.4).  OOO" 
Since  the  dynamics  of  the  critical  states  (5.3)  are  decoupled  from  those  of  the  stable 
states  (5.4),  and  the  dynamics  of  the  stable  states  are  always  LAS  if  ||xc||  is  small  enough, 
the  LASY  of  the  entire  system  (5.3),  (5.4)  (or  equivalently,  (5.1),  (5.2))  is  determined  by 
the  LASY  of  the  subsystem  (5.3).  A  Lyapunov  function  for  the  entire  system  (5.3),  (5.4 
can  be  constructed  if  a  Lyapunov  function  which  determines  the  LAS  of  (5.3)  is  available. 
This  is  summarized  as  follows: 

Theorem  21  When  a  full  normal  form  (5.3),  (5.4)  of  a  system  described  by  (5.1)  and  (5.2) 
exists  and  the  corresponding  transformation  is  given  by  x  =  v(x),  the  LASY  of  the  system 
(5.1),  (5.2),  (or  equivalently,  (5.3),  (5-4)),  is  completely  characterized  by  that  of  (5.3).  If 
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the  subsystem  (5.3)  is  LAS,  and  V c{xc)  is  a  Lyapunov  function  for  the  dynamics  (5.3),  then 
a  Lyapunov  function  for  the  entire  system  (5.3),  (5.4)  is  given  by  V(xc,x3)  =  Vs(x„)  + 
Vc(xc),  where  Vs(xs)  is  constructed  as  in  Lemma  5,  (or  is  any  other  Lyapunov  function  for 
(5.4)).  A  Lyapunov  function  for  (5.1)-(5.2)  is  given  by  V(xc,xs)  =  V(xc,  xs)\x=v(x)  OOO 

Proof:    By  construction,  V(xc,xs)  is  locally  positive  definite.  Consider 

W(xc,xs)  =fV  (xc,xs)  =VC  (xc)+  Vs  (xs)  =f  Wc(xc)  +  Ws(xc,xs). 

By  hypothesis,  in  some  neighborhood  of  xs  =  0;  xc  =  0, 
W(0,0)  =  0 

W(xet0)  =  Wc(xc)  +  Ws(xc,  0)  =  Wc(xc)  <  0;   xc  ?  0 

W(0,xs)  =  Ws(0,xs)<0;  xf#0 

W(xc,xs)  =  Wc(xc)  +  Ws(xc,xs)  <  0;  xc  ?  0,x3  ?  0. 

Therefore,  V{xc,xs)  =  Vs{xs)  +  Vc(xc)  is  a  Lyapunov  function  with  a  locally  negative 
definite  derivative  along  the  dynamics  of  the  entire  system  (5.3),  (5.4),  so  that  (5.3),  (5.4) 
is  LAS. 

Now  taking  V(xc,xs)  =  V(xc,xs)\£=v(x)  as  a  Lyapunov  function  candidate  for  the  sys- 
tem (5.1)  and  (5.2)  yields 

1x7/        \  d'l  dV{xc,xs)     dV(xc,xs)i  w 

W(xc,Xs)  =   —          =   j  \x=v(x)  =  W(xc,Xs)\x=v(x). 


The  local  positive  (or  negative)  definiteness  of  the  function  V(xc,  xs)  (  or  W(xc,  xs))  follows 
since  the  transformation  x  —  v(x)  is  a  local  bijection  (since  the  Jacobian  matrix  at  the  origin 
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is  identity),  hence  the  function  V(xc,xs)  is  a  Lyapunov  function  for  the  system  (5.1)  and 


Although  Theorem  21  gives  a  general  method  for  constructing  a  Lyapunov  function  for 
an  LAS  nonlinear  system  of  the  form  (5.1),  (5.2),  it  is  not  very  practical.  First  of  all,  the 
existence  of  the  full  normal  form  is  not  always  guaranteed.  Secondly,  even  if  a  full  normal 
form  does  exist,  the  derivation  of  this  full  normal  form  requires  one  to  solve  a  partial 
differential  equation.  This  is  impossible  in  general.  However,  for  any  critical  nonlinear 
system,  an  approximate  normal  form  (say  up  to  degree  k)  can  always  be  obtained.  Thus,  it 
is  desirable  to  use  approximate  normal  forms  for  construction  of  Lyapunov  functions.  As  in 
Theorem  21,  the  construction  of  a  Lyapunov  function  for  an  LAS  critical  nonlinear  system 
which  is  in  an  approximate  normal  form  will  depend  on  the  availability  of  a  Lyapunov 
function  for  the  associated  purely  critical  nonlinear  subsystem.  A  general  method  for 
constructing  Lyapunov  functions  for  purely  critical  nonlinear  systems  is  known  only  for 
two  special  cases:  either  the  purely  critical  system  has  one  zero  eigenvalue,  or  it  has  a 
pair  of  purely  imaginary  eigenvalues.  For  LAS  critical  nonlinear  systems  with  such  critical 
dynamics,  Lyapunov  functions  can  be  constructed  by  only  using  approximate  normal  forms 
rather  than  full  normal  forms,  as  will  be  shown  in  next  section. 

5.2    Construction  of  Lvapunov  Functions  for  Nonlinear  Systems  in  Two  Special  Cases 

In  this  section,  the  stability  or  instability  of  nonlinear  systems  in  two  special  cases  is 
studied  using  Lyapunov  functions  and  approximate  normal  forms.  These  two  cases  are: 

CASE  1:  Nonlinear  systems  of  the  form  (5.1)  and  (5.2)  with  Ac  =  0  and  xc  G  R1. 


(5.2) 


000 


/ 


\ 


0  -1 


CASE  2:  Nonlinear  systems  of  the  form  (5.1)  and  (5.2)  with  Ac  = 


and  xc  €  R2. 


\ 


1  0 


/ 
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There  are  many  stability  or  instability  theorems  derived  from  Lyapunov's  theory,  [37]. 
[42].  A  variation  of  the  instability  theorem  which  is  used  intensively  in  rest  of  the  disserta- 
tion will  be  stated  here. 

Theorem  22  [37]  (Instability  Theorem) 

Consider  a  nonlinear  system  x  =  f(x).  Let  x  =  0  be  an  equilibrium  point  for  this  system. 
Let  V  :  D  —>  R  be  a  continuously  differentiable  junction  on  a  neighborhood  D  of  x  =  0, 
such  that  V(0)  =  0  and  V(x)  is  locally  indefinite.  The  equilibrium  point  x  =  0  is  unstable 
ifV(x)  is  locally  definite  (either  positive  or  negative  definite).  000 

5.2.1    CASE  1:  Nonlinear  Systems  with  One  Zero  Eigenvalue 

Consider  a  nonlinear  system  whose  linearization  possesses  only  one  zero  eigenvalue  with 
all  others  being  stable.  An  approximate  normal  form  (up  to  order  k),  where  ak  ^  0  is  the 
first  nonzero  real  number  in  the  sequence  {a2,  a3,  a4,  •  •  •},  has  the  following  structure: 

xc    =    akxkc  +  a(\xc,xs\k+1),  (5.7) 

i,  =  Asxs  +  7s2(^,xs)  +  ---  +  7sfc(^,xs)  +  ^(|xc,xs|fc+1) 

=   Asxs  +  G(xc,xs)xa  +  (3(\xc,xs\k+1)  ,  (5.8) 

where  G(xc,xa)  is  an  ns  x  ns  matrix  function  in  the  components  of  xc  and  xs  with 
G(xc,xs)xs  =  /s2(ic,xa)  +  •  •  •  +  JsixcX,),  and  G(0,0)  =  0.  This  factorization  is  guaran- 
teed since  Jl(xc,xs),  i  =  2,  3,  k  satisfies  7.'(*e»0)  =  0.  The  symbols  a(\xc,x3\h+1) 
and  (3(\xc,xs\k+1)  are  used  to  denote  those  terms  of  degree  k  +  1  or  higher. 

Theorem  23  The  system  (5.7)-(5.8)  is  LAS  if  k  is  odd  and  ak  <  0.  If  ak  #  0  when  k  is 
even,  or  ak  >  0  when  k  is  odd,  the  system  is  unstable.  OOO 
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Proof: 


•  Assume  that  k  is  odd  and  ak  <  0.  Let  Vs(xs)  =  xJPxa,  and  Vc(xc)  =  \x2c,  where 
the  real  matrix  P  =  PT  >  0  is  such  that  (A°)TP  +  PAa  =  Q  <  0.  Let 

F(xc,  x5)  =  Fc(xc)  +  F,(*,)  =  l-x2c  +  xTaPxa. 

The  time  derivative  of  V(xc,xa)  along  the  dynamics  (5.7),  (5.8)  is 

V(xc,xs)    =    akxk+1  +xs{Q  +  H(x)}xs 

+xca(|xc,xs|'c+1)  +  2xJ'P/3(|^,xs|<:+1).  (5.9) 

where  H{x)  d=  G(xc,xs)TP  +  PG(xc,xa).  The  terms  of  degree  k  +  1  or  higher  can 
be  grouped  as  follows: 

xe  a(\xe,  xs \k+i )  +  2  xTs  P  P(\xc,  xs \k+1 )  =  xk+1  *(*„  xs)  +  x*A(xe,xt)x, 

where  $(ic,  i,),  and  A(xc,  xs)  =  AT(xc,xs)  are  some  functions  of  (xe,xa)  with  $(0, 0)  = 
0,  and  A(0,0)  =  0.  The  local  negative  definiteness  of  (5.9)  follows  from  the  equation 
below: 

V  (xc,xa)  =  xk+\ak  +  •(*„*.))  +  xTs{Q  +  H(x)  +  A(i)}  xs. 

•  Assume  that  k  is  odd,  and  ak  >  0.  It  is  clear  that  V(xc,xs)  =  -\x2c  +  xTsPxa  is 
locally  indefinite,  and 

V(xc,xa)    =    -akxk+1  +xs{(Q  +  H(x)}xa 
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-xea(\xc,x,\k+1)  +  2xJ  P  (3(\xc,xs\k+1),  (5.10) 

is  locally  negative  definite.  The  system  is  unstable  when  ak  >  0,  for  k  odd,  by  the 
instability  theorem. 

•  Assume  that  k  is  even  and      ^  0.  It  is  clear  that  V(xc,x3)  =  -akXc  +  xjpxs  is 
locally  indefinite,  and 

V(xc,xs)   =   -a2kxk  +  xs{(Q  +  H(x)}xs 

-ak  a(|xc,i.|*+1)  +  2xTs  P  f3(\xc,xs\k+1),  (5.11) 

with  terms  of  degree  k  +  1  or  higher  be  grouped  as 

-a^a(|xc,x^fc+1)  +  2xfP/3(|xc,Ji|fc+1)  =  ^$(^,xJ  +  xfA(xc,xi)xs, 

for  some  functions  $(zc,xs),  and  A{xc,xs)  =  AT{xc,x3),  with  $(0,0)  =  0,  and 
A(0,0)  =  0.  The  local  negative  definiteness  of  (5.11)  follows  directly  from  this  ex- 
pression, and  the  system  is  hence  unstable  for  k  even. 

OOO 

5.2.2    CASE  2:  Nonlinear  Systems  with  a  Pair  of  Purely  Imaginary  Eigenvalues 

Consider  a  nonlinear  system  whose  linearization  possesses  a  pair  purely  imaginary  with 
all  other  eigenvalues  being  stable.  An  approximate  normal  form  (up  to  order  (2k  +  1)), 
where  m^k+i)  7^  0  is  the  first  nonzero  real  number  in  the  sequence  {7713,  m5,  mj,  •  •  •},  has 
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the  following  structure: 

jcl    =    -uxc2  -  (irf  +  x2c2)n3xc2  +  (x2cl  +  ^2)(i"1)n(2it_i)XC2 

-(x2cl+x2c2)k(n{2k+1)xc2  -  m(2fc+i)Xci)  +  0{\xc,xc2\2k+2)  (5.12) 
icc2    =     uiXd  +  (x2cl  +  x2c2)n3xcl  +  ■■■  +  {x2cl  +  )(*_1)n(2fc-i)*ci 

+(x2cl+x2c2)k(m{2k+1)xc2  +  n{2k+1)xcl)  +  0(\xc,x3\(2k+2)),  (5.13) 
xs    =    ^%  +  L2(xc,x,)  +  ---  +  7i2,:+1)(^,xs)  +  0(|xc,xs|(2fc+2)) 

=    Asx3  +  G(xc,x3)x3  +  0(\xc,xs\2k+2)  ,  (5.14) 

where  m^jt+i)  #  0,  and  G(xc,x3)  is  an  n3  x  n3  matrix  function  in  the  components  of  xs  and 
xc,  with  G{xc,xs)x3  =  J2(xc,xs)  +  ■■■  +  J3(xc,x3),  and  G(0,0)  =  0.  The  representation 
(5.14)  is  due  to  the  fact  that  Jl(xc,0)  =  0,  i  =  2,  3,  •  •  ■,  k.  The  terms  0(\xc,  xs\2k+2) 
denote  terms  of  degree  2k  +  2  or  higher. 

Theorem  21  The  system  (5.12),  (5.13),  (5.14)  is  LAS  if  and  only  if  m2k+i  <  0. 
Proof  of  Theorem  24:  Using  polar  coordinates 

xc2  =  rcosd 
xci    =    r  sin  9 

The  dynamics  (5.12),  (5.13),  (5.14)  is  transformed  to 

r    =    m2k+lr2k+l  +  0{\r,xc2\2k+2) 
%    =    A°x3  +  G(r,x3)x3  +  0(\r,x3\2k+2), 


(5.15) 
(5.16) 
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where  the  dynamics  for  0  have  been  omitted  to  simplify  the  exposition. 

The  rest  of  the  proof  is  similar  to  that  of  Theorem  23,  using  V(r,  xs)  =  ^r2  +  xJPxs 
for  m2k+i  <  0,  and  V(r,xs)  =  —  \r2  +  xJPxs  for  m2k+i  >  0.  When  m2*+i  <  0,  the  system 
(5.12),  (5.13)  and  (5.14)  is  LAS  and  V(xc,xs)  =  \{x2cl  +  x2c2)  +  xjFx,  is  a  Lyapunov 
function.  000 

Example  15  (A  CASE  1  system) 

A  Lyapunov  function  constructed  in  [24]  for  the  study  of  the  LASY  of  the  system 


3 

X\     =     X\X2  +  X1 


i"2    =    —X2  —  2X\ 


is  given  by 


1  7 
V(xux2)  =  x2  +  —x22  +  -x2x2. 


A  normal  form  up  to  degree  3  is  given  by 

i,  =  -x?  +  o(|  -|4) 

i2    =    -x2  +  4x21X2  +  0(\-\i). 
Using  Theorem  23,  a  Lyapunov  function  for  the  system  in  the  normal  form  is  given  by 


V(xc,xs)  =  x\  +  x\. 
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5.3    Lvapunov  Functions  in  the  Original  Coordinates 

Once  a  given  critical  nonlinear  system  whose  critical  dynamics  satisfy  CASE  1  or  CASE 
2  is  transformed  into  an  appropriate  normal  form,  the  stability  (or  instability)  can  be  easily 
determined.  Furthermore,  a  Lyapunov  function  in  the  normal  form  coordinates  (xc,xs) 
can  be  easily  constructed  if  the  system  is  LAS.  This  Lyapunov  function  can  be  used  for 
estimation  of  some  system  performance  characteristics  such  as  stability  region,  convergence 
rate,  and  robustness,  as  well  as  for  other  applications  [24].  However,  it  is  often  preferred 
to  perform  the  stability  and  performance  analysis  on  the  original  system  in  its  original 
coordinates.  A  Lyapunov  function  in  the  original  coordinates  is  hence  desirable. 

Suppose  a  nonlinear  system,  one  of  its  normal  forms,  and  the  associated  transformation 
are  given  respectively  as: 

oo 

x    =    f(x)=  £  fm+\x)/m\  =  Ax  +  f2(x)  +  f3(x)/2\+--- 

m=0 

oo 

j    =    f^)=^fm+1(x)/m\  =  Ax  +  p(x)  +  p(x)/2\+--- 

m=0 

oo 

x    =    u(x)  =  Y,  vm+1(x)/m\  =  x  +  u2(x)  +  u3{x)/2\  +  ---. 

m=0 

From  the  computation  algorithm  of  normal  forms  and  the  associated  transformations  in 
chapter  3,  fk  only  depends  on  /^i)  =  Ax,  f2,  /3,  fk  and  u2,u3,  ■  ■  -,uk.  The  choice 
of  u'+1,t  >  k  does  not  effect  f2,  f3,  and  fk.   To  derive  an  approximate  normal 

form  (say  up  to  degree  k),  a  transformation  can  be  chosen  as  x  =  u(x)  —  x  +  u  (x)  + 

u3(x)/2\  +  h  uk(x)/(k  —  1)!.  To  derive  a  Lyapunov  function  in  the  coordinates  x,  the 

inverse  transformation  x  =  v(x)  -  x  +  v2(x)  +  v3(x)/2\  +  •  •  •  is  needed.  In  general,  the 
inverse  transformation  x  =  v(x)  does  not  have  a  closed  form  even  x  =  u(x)  does.  In  what 
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follows,  we  will  investigate  the  dependency  of  fk  on  =  Ax,  f2,  f3,  •  •  •,  on  fk  and 

v2,  v3,  •  ■  •,  and  we  will  also  derive  an  efficient  algorithm  for  the  computation  of  v2,  v3,  •  •  •  . 

The  homogeneous  function  t'2,  v3,  ■  ■  ■  can  also  be  computed  using  Theorem  13  in  Chapter 
3.  We  first  assume  a  parameterized  inverse  transformation  x  =  v(x)  =  x  +  v2(x)  +  v3(x)/2\  + 
■  ■  ■.  Let 


,m         0°  „m 


n  —  /j 

p(x,a)    =  £pm(*)^=£^+1(*)^, 

m—0  m=0 


m!         '  m! 

m—0  m=0 


Then 


°°  am      00  am 

q(x,a)  =  p(u(x,a),a)  =  £  <lm(x)  — j  =  ^9m+1(i)^, 


m=0  m=0 

is  computed  using  the  algorithm  provided  in  Theorem  13.  By  definition,  q(x,a)  =  x.  Thus, 
a  set  of  equations 

q0(x)  =  x,    q,(x)-0,    for  z  =  1,2,3,--, 

is  obtained.  The  explicit  formulas  for  computing  the  first  few  elements  in  the  sequence 
{v1,  v2,  v3,  v4,  ■  ■  •}  are  given  below: 

q0(x)    =    ux(x)  =  x 

qi(x)    =    0  =  v2(x)  +  U2{x) 

q2(x)    =    0  =  v3(x)  +  ^jf-U2(x)  +  ^f-U2(x)  +  U3(x) 

=  o  =  A-*)  +  d^u2  +  2d-^fu3  +  u4 

d  (  3r_,  ,  dv2(x)TT   ,TT\n,  dqx(x)  dq2(x) 

+¥x  r (x)  +  ~dx~U2  +  U3  U2  +  -dx~U3  +  -dx~U2- 
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Each  vm  for  m  =  1,2, 3,4,  ■••  can  be  derived  by  solving  these  equations.  With  further 
manipulations,  the  solutions  of  vm,s  are  given  below: 

^(x)  =  v0(x)    =  x 

v2{y)  =  v%(y)  =  -Uo(y) 

Ay)  =  Hv)  =  -(^«»)  +  ^W  +  ^(v)) 
«^,)  +  *fW). 

Surprising  as  it  is,  the  computation  of  the  sequence  {vm}  is  simpler  than  that  of  the  sequence 
{um}  when  the  sequence  {U0,  U\,  U2,  •  •  •}  is  known. 
In  general,  we  have 


Recall  that  v0{x)  =  x.  Therefore,  for  a  given  sequence  {Um},  the  sequence  {t;m}  is  uniquely 
determined,  vice-versa.  The  following  theorem  is  obtained  by  observing  equation  (5.17). 


Theorem  25  Let  the  original  system,  one  of  its  normal  forms,  and  the  associated  transfor- 
motion  be  given  respectively  by: 


x   =   f(x)  =  jr  fm+1(x)/m\  =  Ax  +  f2(x)  +  /3(x)/2!  +  •  •  ■  (5.18) 

m=0 
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j    =    /(f)  =  £  fm+1(x)/m\  =  Ax  +  f2(x)  +  p(x)/2\+---  (5.19) 

m=0 

x    =    v(x)=Ytvm+l(x)/m\  =  x  +  v2(x)  +  v3(x)/2\+---.  (5.20) 


m=0 


Then  fk  only  depends  on  /x(z)  =  Ax,  f2,  f3,         and  fk,  and  v2,  v3,  •  •  •,  and  vk.  000 

Finally,  we  have  the  following  theorem  on  the  construction  of  Lyapunov  functions  in 
the  original  coordinates  if  a  given  nonlinear  system  is  LAS. 

Theorem  26  Given  the  system 

x  =  f(x)  =  Ax  +  f2(x)  +  f(x)  +  ■■■,  (5.21) 
its  normal  form,  ( up  to  degree  k ) 

x  =  Ax  +  p(x)  +  ■■■  +  fk(x)  +  ■■■  (5.22) 
is  derived  using  the  transformation 

x  =  v{x)  =  x  +  v2(x)  +  ■■■  +  vk(x).  (5.23) 
IfV{x)  is  a  Lyapunov  function  for  the  LASY  of  (5.22),  then 

V(x)  =  F(:r)U=v(x)  (5-24) 

is  a  Lyapunov  function  for  the  LASY  of  (5.21).  OOO 
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5.4    Examples  and  Discussions 

Several  additional  examples  which  are  calculated  with  the  help  of  Maple  are  presented 
in  this  section.  All  of  the  computations  only  proceed  to  degree  3. 

Example  16  (A  CASE  1  system) 

Consider  the  system: 

2  , 

X1     =     -X1X2  +  xlx3    +  XiX2X3, 

x2  =  -X2  +  x\2  +  xi*32,  (5.25) 
x3    =    -X3  +  x-ix2  . 


A  normal  form  is  given  by 

«i  =  -xl  +  0(\-\% 

x2  =  -x2  +  2xlx2  +  0(\-  |4),  (5.26) 

x3  =  -x3  +  x1x2  +  0(| -|4). 


The  associated  backward  transformation  is 


X1X2X3      xax§      XqX^  xf 

xa  =  Xl+Xlx2-  —  r  +  ~r+j 

~2       ~  ~2   1  =2 


x2    =    x2  +  Xj  -  xix|  +  x\x2,  (5-27) 


=     X3  +  Xi-XiX2. 
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The  associated  forward  transformation  is 


xi3      xij22  ,  XXX2X3  .  XiX32 
xi    =   xi-  x\x2  +  —  H  1  1  — , 

x2  =  x2  -  xi2  +  xi2x2  +  xix32,  (5.28) 
X3    =   x3  -  *i3  +  zi*22  . 


It  is  easy  to  see  that  the  system  (5.26)  by  Theorem  23  is  LAS  using 

K(fi,f3,x3)  =  x\  +  x\  +  xl, 
as  a  Lyapunov  function.  A  Lyapunov  function  for  the  LASY  of  the  system  (5.25)  is  given 

by 


V{x!,x2,x3)    =    \xi-xix2  +  y  +  —2~  +  — 2 —        2~"  j 

,  2 


+  (x2  "  ^l2  +  a'l2;r2  +  XiX32)    +  (*3  -  ^l3  +  *1*22) 


Examvk  17  (A  CASE  2  system) 
Consider  a  Case  2  system: 


ii    =    ~x2  +  xix32  +  XiX2X3  +  Xi2, 


x2    =    xi  +  xa2  +  xxxz2,  (5.29) 


X3    -     -x3  +  xax2  . 


A  normal  form  up  to  degree  3  is  given  by 


_•  _      (Sx1-7x2)(xl  +  xl)  4 

x\  =  -x2  —  +  U(\-\), 

*2  =  xi  ^2  +  u(\-\)' 

i3  =  -^3  +  0(1  •  I4)  ■ 


The  associated  backward  transformation  is  given  by 


x\     2x\x2     2x\     1  x\     x\x2  x\xz 
xi    =    *l-  z+—$  3~  +  l8        2  5~~ 

29£i£|  _  3  x-ix2x3  _  xxx\  _  11  x3  ,  *2*3  x2x\ 
+    36  10  2         12         5  4 

x\     2xxx2     2x\     x\     \lx\x2  x\x3 
x2    =    x2  +  J  +  —  +  —  +  l2  +  ~36~  +  lT 

_  x\x2X3  _  X\x\  _  x?,x3     2_i|  2i!^ 
10  4         10        9  3 

ii  ,  x\x~2  .x\     2x\     x\x2  ,  4x1x2!  2x\ 


The  associated  forward  transformation  is  given  by 


xi2     2xix2     2x2      xi3     X\2Xj  t  X\2xz 

xi  =  ^  +  T'~T~  +  ^~  +  l8~~T'  +  ~r 

13 xix22     3x1X2X3  ^  X!X32     5  x23     x22x3  x2x3 


36  10  2          12  5 

x\2  2xix2  2x2      xi3     29xx2x2  X!2x3 

x2    =    x2-—z  3  3~~_l2"         36      ~  10 

2xix22  X1X2X3  X1X32     2x23  x22x3 

+     3  10  4910' 
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The  system  (5.30)  is  LAS  by  Theorem  24.  A  Lyapunov  function  for  (5.30)  can  be  chosen 
as 

V(xi,x2,  x3)  =  x\  +  x\  +  x\  . 
A  Lyapunov  function  for  the  LASY  of  the  system  (5.29)  is  given  by 

/        xi2     2xyx2  |  2x22  ,  23  x:3     W12  xrt 
V{xux2fx3)    -    yx1  +  -  —  +    3    +    36  12     +  5 

2xix22     3xix2x3     xqx32     2x23  _  x22x3  x2x32\2 
+      9      +      10  2~3~5  4  J 

,  /         X!2      2xax2      2  x22      XX3      25Xq2X2  Xi2X3 

+  (I2_t-^  r  +  i"+^8  ir  (5-33) 

llXiX22      XiX2X3  +  XqX32  +  29  x23  +  X22" 


— V 


12  10  4  36  10 

/        Xi2     xax2     x22  xi3\2 

+  {X3+  —  --r--^--5-)  ■ 

The  expression  for  F(xi,x2,i3)  can  also  be  written  out  explicitly.  This  expression  is 
intentionally  omitted  due  to  its  complexity.  The  local  negative  definiteness  of  Vr(x1,x2,x3) 
is  guaranteed  by  Theorem  24.  Since  the  Hessian  of  F(x!,x2,x3)  at  the  origin  is  negative 
semi-definite,  the  local  negative  definiteness  of  V(xux2,  x3)  is  quite  difficult  to  verify  using 
other  methods. 


Example  18  (A  CASE  2  system) 

Consider  a  simple  two-dimensional  Case  2  system: 


xi    =    -x2  +  x2 

x2    =    xi  +  x2  .  (5.34) 
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A  normal  form  up  to  degree  3  is 


xi    =  -*2-f  +  +       +  0(|-|4) 

=  -x2-(i^lMM±fI)+0(|.|4) 

*2     =  JH-^-fl^-lM.M  +  0(|.|4) 

=  ,i_(3x2  +  7x1Hxf  +  xl)+Q([  |4)  (5  35) 


The  associated  backward  form  transformation  is  given  by 


_x\     2x;x2     2xj     23  x^     3x1     2x\  x\x2 
Xl_Xl~T+3~3  36  4        9  3 


.     ,   Xj   ,  2X!X2   ,  2jj       X?    ,   5X^X2       *1*2    i   *2  ,,0^ 


The  associated  forward  transformation  is  given  by 


xa2     2xix2     2x22     2xi3     xi2x2     7xix22  7x23 


Xl    ~    Tl+  3         3     +    3    +    9  3  36  12 

_  fi!  _  2x1x2  _2^V     35  xt2x2  jV       (g  3?) 

T2    -   12      3         3  3    +  12  +      36     +      6     +    18        1  ' 


The  LASY  of  the  system  (5.37)  is  evident  by  Theorem  24.  A  Lyapunov  function  is 
given  by 

V(xi,x2)  =  x\  +  x\  . 


A  Lyapunov  function  for  the  system  in  the  original  coordinates,  (5.34),  is  given  by 


(        xj2     2xtx2     2x2*  ,  2xic 
V(xa,x2)    =      x1  +  1  _  +  _  +  _ 


Xi2X2 


7xix2" 
36 


12 


7^ 


/  x-i2     2xix2     2x2  xi3 

[+X2-~3 — 3  r  + it 

t  35x^X2  ,  5xxx22  ,  7x23V 

+  ^6_  +  _6_  +  lT      •  (5-38) 


The  local  positive  definiteness  of  V(x\,x2)  is  obvious. 

The  time  derivative  of  V(xi,x2)  along  the  dynamics  (5.34)  is  given  by 


dV         247xx6x2     173  Ji6     19  xt7     7x26     10  x25 
~dt  108  216         54         8  9 

32xi3x22     61  xx2x23     413x25X!2  _  212j13j23 
+      9       +       9       +       108  27 

983xi5x22     1703xi2x24     1315x!4x23  29xix2 
+      108  216       +       108  36 

439xx3x24  _  641x!5x2  _  661  xt4 j22  8x^ 
+      54  108  72     "  9 


5 


+2x14*2-^-^-x12x22.  (5.39) 


Notice  that 


^  =  -^I  +  ^)2  +  0(|-|5)-  (5-40) 


In  this  case,  the  local  negative  definiteness  of  (5.40)  is  also  evident  by  using  polar  coordi- 
nates. 

This  example  has  been  studied  by  Mees  and  Chua  in  [44],  as  well  as  by  Fu  and  Abed 
in  [24].  In  [24],  a  Lyapunov  function  for  (5.34)  is  given  as 


V-(*x,x2)  =  \{x\  +  x2)  -  x\x2  +  d_M  +  ^_*±l_*A.  (5.41) 
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The  time  derivative  of  V(xi,x2)  along  the  dynamics  (5.34)  is  given  by 


^  =  -\(*l  +  4)2  +  0(\-\5).  (5.42) 


It  appears  that  (5.41)  is  simpler  than  (5.38)  in  some  sense.  In  what  follows,  it  will  be  shown 
that  this  difference  is  only  superficial.  The  method  proposed  in  this  paper  can  generate  a 
Lyapunov  function  which  is  identical  to  (5.41),  as  illustrated  below. 

As  we  mentioned  before,  the  normal  form  for  a  given  system  is  not  unique.  It  is  easily 
verified  that  the  approximate  normal  form  (5.35)  can  also  be  derived  using  the  following 
transformation 


Xl2     2xix2  ,  2z22  ,  (l/9-«3,2,4)Ji3 

XI         =         *X  +   -g  —    +—   +  2 

(-l/3  +  a3,2,3)*l2S2    ,    (-13/18-03,2,4).  2 

+  2  +  2  12 


,  (-5/6  +  03,2.3)  *23 
+  2 


xi2     2xxx2     2x22  ,  (-1/6  -  03,2,3) xic 
i2    =    z2  r  V 


3         3  3  2 

(29/18-o3,2,4)  2  ,  (4/3- 03,2,3)  xix22 
+  ~2  *2  +  ~2  

+  (4/9  -03,2,4)^  j  (5  43) 

2 


where  03,2,4  and  03,2,3  are  arbitrary  real  numbers.  The  forward  transformation  (5.37)  is 
obtained  by  setting  a3,2,4  =  03,2,3  =  —1/3.  Choosing 


V(xl,x2)  =  -  (x\  +  xty  , 


so 


a  parameterized  Lyapunov  function  Vp(xi,x2)  for  the  system  (5.34)  derived  using  the  for- 
ward transformation  (5.43)  is  obtained  by  the  substitution.  Notice  that  no  matter  what 
values  03,2,4  and  03,2,3  take, 

^  =  -\(xl  +  4)2  +  0(\-\*).  (5.44) 

The  degree  of  the  lowest  order  homogeneous  polynomial  in  (5.44)  is  four,  and  (5.44)  is  locally 
negative  definite  if  its  lowest  order  homogeneous  polynomial  is  negative  definite.  Therefore, 
a  modified  parameterized  Lyapunov  candidate  Vnew(xi,  x2),  obtained  from  Vp(xi,x2)  by 
truncating  terms  of  degree  5  or  higher,  is  as  follows, 

l/0      ,\        9        2i23     x\3  xi3x2 
Vne"(xux2)   =    -(xl  +  xj)-x12x2--^-  +  -±  Y- 

(-2a3l2,4  +  8/3)x12a:22  xxx2z 
+  2  4 

(1/3- 03.2,4)  *i4     (-03,2.4  + 4/3)  x24 
+  2  2 

The  local  positive  definiteness  of  Vnew(xi,  x2)  is  also  evident.  The  Lyapunov  function  (5.41) 
is  obtained  by  setting  a3,2,4  =  4/3  in  Vnew.  The  value  of  a3,2,3  does  not  matter. 

Remark  9  Actually,  all  Lyapunov  functions  constructed  in  the  above  examples  can  be  sim- 
plified. This  topic  will  be  pursued  further  in  the  next  chapter.  OOO 

5.5  Summary 

In  this  chapter,  we  first  presented  a  general  algorithm  for  construction  of  Lyapunov 
functions  for  nonlinear  critical  systems  when  they  are  LAS.  It  was  shown  that  when  a  given 
critical  nonlinear  system  is  transformed  into  a  full  normal  form  (assume  it  is  existed),  a 
Lyapunov  function  of  the  entire  system  in  the  normal  form  can  be  chosen  to  be  the  sum  of 
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the  individual  Lyapunov  functions  corresponding  to  the  critical  and  the  stable  dynamics, 
respectively.  A  Lyapunov  function  for  the  system  in  the  original  coordinates  can  be  derived 
by  substitution  of  the  appropriate  transformation.  In  the  Section  5.2  and  5.3,  an  explicit 
algorithm  for  the  construction  of  Lyapunov  functions  for  two  classes  of  critical  nonlinear 
systems  was  provided.  This  construction  procedure  only  required  an  approximate  normal 
form  up  to  a  certain  degree.  The  Section  5.4  provided  several  examples  to  highlight  use  of 
the  algorithm.  The  issue  of  simplification  of  Lyapunov  functions  is  mentioned,  and  will  be 
pursued  further  in  Chapter  6. 


CHAPTER  6 
SIMPLIFICATION  OF  LYAPUNOV  FUNCTIONS 

In  Chapter  5,  we  proposed  an  explicit  algorithm  for  the  construction  of  Lyapunov 
functions  of  nonlinear  systems  in  either  CASE  1  or  CASE  2.  In  this  chapter,  we  will 
provide  a  method  to  refine  those  constructed  Lyapunov  functions.  As  we  mentioned  before, 
the  construction  of  Lyapunov  functions  for  nonlinear  systems  in  either  CASE  1  or  CASE 
2  has  been  studied  by  some  other  researchers.  A  comparison  among  Lyapunov  functions 
constructed  using  different  methods  motivated  the  work  presented  in  this  Chapter.  Some 
heuristic  examples  are  given  in  Section  6.1.  In  Section  6.2,  we  present  the  method  for 
refinement  of  Lyapunov  functions.  Several  examples  are  given  in  Section  6.3  to  illustrate 
use  of  the  method.  This  chapter  is  concluded  by  a  brief  summary  given  in  Section  6.4. 

6.1    Heuristic  Examples 

Example  19  (A  CASE  1  system) 
The  system 

ii    =    xix2  +  xl 

x2    =    -x2-2x\  (6.1) 
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is  considered  by  Fu  and  Abed  in  [24].  The  origin  is  an  LAS  equilibrium  point  by  using 
center  manifold  theory.  A  Lyapunov  function  constructed  in  [24]  is 

1  7 
V(xi,x2)  =  x\  +  —x]  +  -x\x2. 

OOO 

Furthermore,  Fu  and  Abed  proved  that  if  a  CASE  1  system  is  LAS  and  a3  <  0  is  the 
first  nonzero  real  number  is  the  sequence  {a2,  03,  a4,  •  •  •},  a  Lyapunov  function  including 
only  quadratic  and  cubic  terms  can  be  constructed.  Similarly,  if  a  CASE  2  system  is  LAS 
and  m3  <  0  is  the  first  nonzero  real  number  is  the  sequence  {7713,  7715,  m.7,  ■  •  •},  a  Lyapunov 
function  including  only  quadratic,  cubic  and  quartic  terms  can  be  constructed.  The  proof 
of  their  method  is  based  on  the  following  proposition. 

Proposition  1  [24]  (Local  Definiteness  of  a  Class  of  Bivariate  Functions):  The  scalar  func- 
tion 

A(x,  y)   =   a2oX2  +  a30a;3  +  a21  x2y  +  ai2  xy2  +  a40x4 

+a3ix3y  +  a22  x2y2  +  ai3xy3  +  a04y4  +  0(\x,  y\5) 

in  the  real  variables  x  and  y  is  locally  negative  definite  in  an  open  neighborhood  of  the  origin 
if  a2o  <  0,  a04  <  0,  and  a\2  -  4a20a04  <  0.  OOO 

If  determining  the  LASY  of  a  nonlinear  system  needs  the  consideration  of  nonlinear 
terms  of  degree  4  or  higher,  Fu  and  Abed  did  not  give  any  clue  for  the  construction  of 
Lyapunov  functions.  On  the  contrast,  the  algorithm  we  proposed  in  chapter  5  is  applicable 
to  any  CASE  1  or  CASE  2  nonlinear  system  no  matter  how  complicated  the  nonlinear 
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system  is.  If  a  nonlinear  system  is  LAS  and  a-ik+i  <  0  is  the  first  nonzero  real  number  in  the 
sequence  {a2,  a3,  a4,  •  •  •}  for  CASE  1,  (or  m2fc+i  <  0  is  the  first  nonzero  real  number  in  the 
sequence  {m3,  m5,  m7,  •  •  •  }  for  CASE  2),  Lyapunov  functions  in  the  original  coordinates 
(xc,  xs)  which  include  terms  of  degrees  2  to  2(2k  +  1)  can  be  constructed  constructed  using 
the  algorithm  provided  in  Chapter  5  The  idea  of  constructing  a  simpler  Lyapunov  function 
for  a  given  nonlinear  system  inspired  the  research  presented  in  this  chapter.  We  want  to 
prove  that  it  is  possible  to  construct  Lyapunov  functions  in  the  original  coordinates  (xc,  xs) 
that  only  include  terms  of  degree  from  2  to  2k  +  2  ( instead  of  2(2k  +  1) ). 

Examvle  20  (A  CASE  1  system) 

In  Chapter  5,  a  Lyapunov  function  has  been  constructed  for  the  following  CASE  1 
system  described  by 

x    =    -xy  +  xz2  +  xyz, 

y  =  -y  +  x2  +  xz2,  (6.2) 
z   =    -2  +  xy. 

Choosing  a  Lyapunov  function  for  the  system  in  the  normal  form  as 

V(x,y,  z)  =  x2  +  y2  +  z2, 

a  Lyapunov  function  for  the  system  in  the  original  coordinates,  constructed  using  the 
method  proposed  in  Chapter  5,  is 

...         .  (  x3     xy2     xyz  xz2V 

V(x,y,z)   =    [*-*y+Y  +  -T  +  ^-  +  —) 
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+  (y  -  x2  +  x*y  +  x,2)2  +  (.  -  x3  +  xy2)2 .  (6.3) 
The  time  derivative  along  the  trajectory, 

—  =  -2  z2  -  2 1/2  +  4  x2y  +  2  xt/2  +  2  x32  -  4  x4  -  6  xy2z  -  4  xyz2  +  0(\  ■  |5)  (6.4) 
dt 

is  guaranteed  to  be  negative  definite  in  at  least  a  small  open  neighborhood  of  the  origin. 
Due  the  complexity  of  dV/dt  here,  the  terms  of  degree  5  or  higher  are  omitted.  In  the 
following  we  will  prove  that  a  function  obtained  by  truncating  the  terms  of  degree  5  or 
higher  from  V(x,y,  z)  in  (6.3)  is  still  a  Lyapunov  function  for  the  system  (6.2).  First,  we 
need  to  prove  the  following  Lemma. 

Lemma  6  The  function  A(x,y,z)  =  -2z2-2y2+4x2y+2xyz+2x3z-4x4-6xy2z-4xyz2 
is  locally  negative  definite.  Any  other  analytic  function  A{x,y,z)  generated  by  A(x,y,z) 
such  as  A(x,  y,  z)  =  A(x,  y,  z)  +  0(\  ■  |5)  is  also  locally  negative  definite.  OOO 

Proof:  The  function  A(x,y,z)  can  be  re-written  as  follows: 

A(x,y,z)   =    -2  z2  -  2  y2  +  4  x2y  +  2  xyz  +  2  x3z  -  4  x4  -  6  xy2z  -  4  xyz2 


=    —  a\y2  +  2xyz  —  a^z2  (6-5) 

-a3y2  +  4x2y  -  a4x4  (6.6) 

-z2(a5  +  4xy)  (6.7) 

-y2(a6  +  6xz)  (6.8) 

-a722  +  2x32  -  a8x6  (6.9) 

-x4(a9-a8x2)  (6.10) 
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where  ai  -  a9  are  some  positive  real  numbers.  Recall  that  ax2  +  2bxy  +  cy2  is  negative 
definite  if  and  only  if  c  <  0,  b  <  0  and  b2  -  ac  <  0.  Therefore,  if  the  conditions 

1.  x2  -  aia2  <  0, 

2.  4  -  a3a4  <  0, 

3.  6a  =  a5  +  4xy  >  0, 

4.  62  =  a6  +  612  >  0, 

5.  1  -  a7a8  <  0, 

6.  63  =  a9  -  a8i2  >  0 

hold  in  a  small  neighborhood  of  origin  (1  =  y  =  z  =  0)  subject  to 

1.  ai  +  a3  +  a6-  2, 

2.  a2  +  a5  +  a7  =  2. 

3.  a4  +  a9  =  4, 

A(z,y,  2)  will  be  locally  negative  definite.  These  conditions  are  absolutely  achievable  pro- 
viding x,  y,  and  z  are  small  enough.  For  example,  all  inequalities  hold  if  we  can  choose 


1  ,  3  7  1  1  1  A  1\ 

(ai  =  -,a2  =  l,a3  =  -,a4  =  -,a5  =  -,a6  =  -,a7  =  -,a8  =  4,a9  =  -) 


and       |y|,  \z\  <  ^.  The  higher  order  terms,  0(\  ■  |5),  can  be  re-written  as 


0(|  •  |5)  =  y2hi(x,  y,  x)  +  yzh2(x,  y,  z)  +  z2h3(x,  y,  z)  +  x4h4(x,  y,  z). 
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where  hi  -  h4  are  some  polynomials  vanishing  at  origin.  Consequently,  we  can  re-write 
A(x,y,z)  as 

A(x,y,z)   =    {-aiy2  +  2xyz  -  a2z2) 
+(-a3j/2  +  4i2y  -  aAxA) 
+(-a7z2  +  2x3z  -  a8x6) 
-z2  {{a5  +  4xy)  -  h3(x,  y,  z)}  +  yz  h2(x,  y,  z) 

—y2  {(ae  +  612)  -  hi(x,  y,  2)} 
-x*{(a9  +  asx2)  -  h4(x,  y,  z)}. 

The  function  A(x,  y,z)  is  indeed  locally  negative  definite.  OOO 
With  this  lemma  a  new  Lyapunov  function,  Vnew,  can  be  derived  from  (6.3)  by  trun- 
cating the  terms  of  degree  5  or  higher.  We  have 

Vnew(x,y,z)   =    x2 +  y2  +  z2 -4x2y  +  2x4 -2x3z  + z2x2 +  4x2y2 

+yzx2  +  2  xyz2  +  2  xy2z.  (6.11) 

One  drawback  of  this  approach  is  that  this  Vnew(x,  y,  z)  is  only  locally  positive  definite, 
whereas  V(x,  y,  z)  in  (6.3)  is  positive  semi-definite  globally. 

Later  we  will  show  that  a  Lyapunov  function  for  (6.1)  which  only  includes  quadratic 
and  cubic  terms  can  also  be  constructed  using  the  method  proposed  in  the  next  section. 
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6.2    Main  Results:  Simplification  of  Lvapunov  Functions 

Due  to  the  similarity  between  CASE  1  and  CASE  2  systems,  only  CASE  1  systems  are 
considered  in  this  section.  Consider  a  given  nonlinear  system  in  CASE  1, 

xc  =  Acxe  +  fe(xc,xt),  xceRnc,  (6.12) 
x,    =    A3x3  +  f3(xc,x3),   x3eRn>.  (6.13) 


The  LASY  of  the  entire  system  is  decided  by  a  set  of  real  numbers  {a2, 03,  a4,  as,  •  •  •}.  This 
was  discussed  in  Chapter  5.  If  such  a  critical  system  is  LAS,  then  the  first  nonzero  number 
in  sequence  {a2,  a3,  a4,  as,  •  •  •}  must  be  negative  with  odd  subscript,  say  a2k+i  <  0.  The 
normal  form  up  to  degree  2k  +  1  is  given  by 

Zc    =  a2k+1x2ck+1  +  o{\xc,x3\2k+2)  (6.14) 

x3    =  A  is -(- /3(xc,  xs) 

=  Asx3  +  fl(xe,  x.)  +  •  •  •  +  7f +2(ic,  i.)  +  0(\xct  x3\2k+2) 

=  Asx3  +  e(xc,x3)x3  +  0(\xc,xs\2k+2)  (6.15) 


 ■  — 2   _  _ 

where  fs(xc,x3  =  0)  =  0,  *  =  1,  2,  •  •  • ,  2k  +  1,  and  0(xc,xs)x3  =  f3(xc,xs)  +  ---  + 

/,  (Xc,X3). 

Let  P  =  PT  >  0  be  a  matrix  such  that  AsTP  +  PA"  <  0.  By  Theorem  23,  a  Lyapunov 
function  for  the  system  (6.14)  and  (6.15)  can  be  chosen  as 


V(xe,x,)=  -x2e  +  xJPxa  (6.16) 
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Let  the  forward  transformation  be 

x  =  v(x)  =  x  +  v2(x)  +  ■■■  +  v2k+1(x).  (6.17) 

Using  (6.16),  a  Lyapunov  function  for  the  system  (6.12)  and  (6.13)  (in  the  original  coordi- 
nates (xc,xs) )  has  the  form: 

V(xc,  xs)    =    V  (x)  \x=v(x) 

=    V2(xc,  x,)  +  •  •  •  +  V2k+2(xc,  x.)  +  •  •  •  +  V4k+2(xc,  xs).  (6.18) 

Let  the  function  Vnew(xc,xs)  be  obtained  from  V(xc,xs)  by  cutting  terms  V2k+3(xc,xs)  + 
...  +  V4k+2(xc,xs),le., 

Vnew(xc,  xs)  =  V2(xc,  x.)  +  •  •  •  +  F2fc+1(xc,  xs)  +  V2k+2(xc,  xs).  (6.19) 

We  want  to  prove  the  following  theorem. 
Theorem  27  The  function 

Vn™(xc,  x.)  =  V2(xc,  x.)  +  •  •  •  +  V2k+\xc,  x.)  +  V2k+2(xc,  xs)  (6.20) 

is  a  Lyapunov  function  for  the  system  (6.12),  (6.13).  OOO 

In  order  to  prove  this  theorem,  we  need  to  introduce  a  key  lemma  first. 

Consider  a  given  function  which  is  given  as  a  summation  of  homogeneous  polynomials: 


n(z)  =  n2(x)  +  •  •  •  +  nk(x),  xeRn. 


(6.21) 
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Assume  that  II(i)  is  locally  positive  definite,  and  furthermore,  any  analytic  function  II(x) 
of  form 

E(x)  =  n2(x)  +  •  •  •  +  Uk(x)  +  0(\x\k+1),  (6.22) 

is  also  locally  positive  definite,  where  0(|i|fc+1)  denotes  any  combinations  of  terms  of  degree 
k  +  1  or  higher.  Let 

x  =  v(y)  =  y  +  v2(y)  +  v3(y)  +  ■■■ 
be  an  analytic  transformation.  It  is  obvious  that 

A(y)    ='   n(x)|x=v(y)  =  A2(j/)  +  .."•  +  A*(y) +A*+1(y)  +  ■  •  • 

=}   A(»)  +  a(||y||fe+1)  (6.23) 

is  still  analytic  and  locally  positive  definite.  Here,  a(|t/|fc+1)  denotes  some  combinations  of 
terms  of  degree  k  +  1  or  higher.  We  claim  that  the  following  lemma  holds. 

Lemma  7  The  function  A(y)  =  A2(y)  +  ■■■  +  Ak(y)  derived  from  (6.23)  is  locally  positive 
definite.  Furthermore,  any  analytic  function  A(y)  =  A(y)  +  0(\y\h+1)  =  A2(y)  +  ■■•  + 
Ak(y)  +  0(\y\k+1)  is  also  locally  positive  definite.  OOO 

Proof:  Let  the  exact  inverse  of  the  analytic  transformation  x  =  v(y)  =  y+v2(y)  +  v3(y)-\  

be 

y  =  u(x)  =  x  +  u2(x)  +  u3(x)  +  ■  ■  ■ . 
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By  the  inverse  function  theorem  [29],  the  analyticity  of  v(y)  at  the  origin  y  =  0  guarantees 
that  of  u(x)  at  the  origin  x  =  0.  From  II(i)  |x=v(y)=  A(y)  +  a(\y\k+1),  we  have 

A(y)  =  n(x)  -  a(\y\k+1)  =  U(x)  -  /3(|z|fc+1)  (6.24) 

for  some  function  /?(•).  Thus  the  local  positive  definiteness  of  function  A(y)  follows  from 
the  local  definiteness  of  function  n(x). 
Notice  that 

AG,)   =   A(y)  +  0(\y\k+1)  =  Il(x)  +  0(\y\k+1)  -  P(\x\k+1) 

=    n(x)  +  7(|x|fc+1)-/?(k|fc+1)  (6.25) 

for  some  function  7(-).  The  local  positive  definiteness  of  function  A(y)  follows  from  the 
local  definiteness  of  function  II(z)  consequently.  OOO 

We  are  now  ready  to  prove  the  Theorem  27. 
Proof  of  the  Theorem  27:  Differentiating  both  sides  of  (6.18)  with  respect  to  time  along  the 
trajectory  of  the  CASE  1  system  (6.14)  and  (6.15)  yields 

V{xc,xs)    =  V(xc,xs)\x=v(x) 

=    [a2k+i  x2ck+2  +  xTs  [{A-  +  Q(xc,  xs))TP  +  P(A~  +  Q(xc,  z.)j)  xs 


W2(xc,  *.)  +  .-+  W2k+2(xc,  xs)  +  0  (l^|2fc+3)]_=v 
=   W2(xc,  x.)  +  •  •  •  +  W2k+2(xc,  xs)  +  O  (\x\2k+2)  . 
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We  have  shown  the  local  definiteness  of  both  W  (xc,xs)  +  •  •  •  +  W  {xc,xs)  and 
W2{xc,xs)  +  ---  +  W2k+2(xc,xs)  +  0  (|x|2fc+3)  ■  The  local  definiteness  of  both  W2(xc,xs)  + 
■  ■  ■  +  W2k+2(xc,  x.)  and  W2(xc,  xs)  +  --  ■  +  W2k+2(xc,  x3)  +  0  (\x\2k+3)  Mows  directly  from 
Lemma  7. 

Since  the  terms  V2k+3(xc,xs)  +  h  V4k+2(xc,x3)  in  V(xc,xs)  have  no  influence  on 

the  local  negative  definiteness  of  V(xc,xs),  these  terms  can  be  truncated.  Of  course, 
Vnew(xc,  xa)  is  still  locally  positive  definite  since  its  Hessian  is  a  positive  definite  matrix. 

OOO 

6.3  Examples 
Example  21  (Continuation  of  example  19) 

The  following  normal  form  transformation  for  the  system  (6.1) 

Xi    -    X\  +  c\x\  -  X\X2  +  cAx\  +  (-ci  -  b2)  x\x2  +  c&x\i\, 
x2    =   x2  -  2x\  +  b2xix2  +  (1/2  -  c6)x2 
-4caz?  +  b5x\x2  +  Mi*2  +  Ma, 

yields  a  normal  form  up  to  degree  3 

*2   =   -*2  +  (1/2  -  ce)xl  +  (66  -  b2(l/2  -  c$)xxx\ 
+4  x\x2  +  (2  b7  -  2  (1/2  -  c6)2)  if, 
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where  the  coefficients  {ci,  C4,  ce,  62^  h,  be,  h}  are  free  parameters.  The  forward  transforma- 
tions is 

2 

9  /  n  \         o  9  J  1  J  1 

xa    =    xx  -  c^j  +  xix2  +  [2  +  2ci  —  C4J  ij  —  2cjXjX2  4  2  1 

*2    =    *2  +  2xx2  -  62^i^2  +  (^6  -  1/2) x22  -  262^i3 
+  (2  +  622  +  ca62  +  4  c6  -  65)  Xi2I2 

+  (-3  c662  +  -j-&ej  *i*22  +  (-2  c6  +  1/2  +  2  c62  -  fc7)  x23. 
Choosing  a  Lyapunov  function  for  this  system  in  normal  form  as 


x2 

V(xi,x2)  =  x\  +  -|, 


gives  a  Lyapunov  function  in  the  (xi,i2)  coordinates  as,  after  truncating  terms  of  degree 
5  or  higher, 


Vnew,       ,x     _     T2      *2      ?  .  _3  ,   9  *1*2       *2fj*l       /  C6  1 

3 


+  ( j  +  5  Cl2  -  2  c4)  *?  +  (-6  Cl  -  ^  x2x? 


+  l~8~  +   16  +    8    +  4  ~  8  r2  1 


.2„2 


3c662  .  h     b6     b2  (c6  -  l/2)\  3 
+  ^"l~  +  l6-  J  8  )X2Xl 


(  c6      1      c62     67     (c6-l/2)2\  4 

+  rT  +  ^  +  T_¥+    16  r2> 


A  very  simple  Lyapunov  function 


F-(x1,x2)  =  x2  +  g  +  ^2 
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is  generated  by  setting  {cj  =  0,c4  =  17/8,  c7  =  1/2,  b2  =  0,b5  =  20, 66  =  0, 67  =  —1/2}. 
The  local  negative  definiteness  of 

dVnew(xi,x2)  _     x2y     5x4     9x2y2     9x4y  y2 
dt  2   ~    2   +     2     +    2    ~  8  ' 

can  be  easily  verified  Proposition  1.  OOO 
Example  22  (Continuation  of  Example  17  in  Chapter  5) 

Using  Theorem  27,  a  simplified  Lyapunov  function  for  the  CASE  2  nonlinear  system  (5.29) 
is  given  by: 

o       0       9     2xi3          ,        2xi2x3  2xix2x3 
Vne-(x1,x2,x3)    =    x12  +  x22  +  x32  +  -j--2x12x2  +  ^-^  

2x22x3     4x23     77  xx4     29x2X!3     124xi2x22  2x3Xj2x2 
5  3~ +    50    ~     50     +       25       +  5 

,    o     29xjx23     x3xix22     127  x24  x3x23 

+    13  ii  + — —  =  +  — TR — +  ~ T~  ■ 

50  5  50  5 

It  is  also  easy  to  verify  the  local  negative  definiteness  of 


dv-idt  =  -2x32  +  i^  +  i^-i^  +  i^i! 

0              5              0  25 
4x!X23     29  x24     23x12x22     4X3X13     29xx4  5 
 25  50  2T-  +  ^  50~  +  O(W  }- 


6.4  Summary 

In  this  chapter,  the  procedure  for  the  construction  of  Lyapunov  functions  for  nonlinear 

systems  was  refined.  If  the  local  stability  of  a  given  nonlinear  system  x  =  Ax  +  /2(x)H  h 

f2k+1(x)  +  •  •  •  is  equivalent  to  that  of  the  system  x  =  Ax  +  /2(x)  +  1-  /2*+1(x),  it  was 
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shown  that  polynomial  Lyapunov  functions  which  include  only  terms  of  degree  2  to  (2k+2) 
can  be  constructed.  Since  normal  form  transformations  are  not  unique,  the  constructed 
Lyapunov  functions  can  be  further  simplified  by  carefully  choosing  those  free  parameters  in 
the  normal  form  transformations.  Several  examples  were  provided  to  illustrate  use  of  the 
refined  algorithm. 


CHAPTER  7 

CONSTRUCTION  OF  LYAPUNOV  FUNCTIONS  -  II 


In  this  Chapter,  we  will  provide  an  alternative  algorithm  for  the  construction  of  Lya- 
punov  functions  for  critical  CASE  1  and  CASE  2  critical  nonlinear  systems.  This  algorithm 
requires  that  only  the  dynamics  of  critical  states  be  transformed  into  a  normal  form  up  to 
an  appropriate  degree.  In  Section  7.1  a  heuristic  example  is  given.  This  example  partially 
motivated  the  algorithm  proposed  in  this  chapter.  Section  7.2  provides  the  main  theorem 
of  this  chapter:  a  method  for  the  construction  of  Lyapunov  functions  for  any  CASE  1  or 
CASE  2  nonlinear  system  once  it  is  in  an  appropriate  partial  normal  form.  In  Section  7.3, 
Lyapunov  function  in  the  original  coordinates  are  derived,  and  simplified.  The  computation 
of  partial  normal  forms  is  discussed  in  Section  7.4.  More  examples  are  presented  in  7.5  for 
illustration  of  the  algorithm.  A  brief  summary  in  Section  7.5  concludes  this  chapter. 


7.1    A  Heuristic  Example 


Example  23  (A  heuristic  example) 


Consider  the  system 


3  ,     2  2 

x    =  —x+xy 


(7.1) 


y    =    -y  +  x2  +  y3. 


(7.2) 
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Since  any  center  manifold  must  have  the  form  y  =  h(x)  =  0(\ •  |2),  we  can  conclude,  without 
actually  calculating  the  center  manifold,  that  the  dynamics  along  the  center  manifold  is 

i  =  -e+o(\-\3) 

which  is  locally  asymptotically  stable.  Therefore,  the  entire  system  (7.1)-(7.2)  is  LAS. 
Up  to  degree  3,  an  approximate  normal  form  of  (7.1)-(7.2)  is  given  by: 

x   =   -x3  +  0(H4) 

y  =  -y  +  0(\-\4) 

with  the  associated  transformations  given  by 


x    =    x  +  0(H4), 

y  =  y  +  x2-j  +  0(\-\% 


and 


x    =   x  +  0(|  •  |4) 

y     =     y-Xl+yl  +  0(\.\<). 


The  subsystem  (7.1)  can  be  viewed  as  in  an  approximate  normal  form,  but  the  subsystem 
(7.2)  is  not.  Nevertheless,  without  further  applying  either  center  manifold  or  normal  form 
theory,  the  LASY  of  (7.1)-(7.2)  is  evident.  This  observation  provides  an  inspiration  for 
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the  major  results  -  the  construction  of  Lyapunov  functions  using  partial  normal  forms, 
presented  in  this  chapter. 

7.2    Main  Results 

For  simplicity,  consider  an  analytic  nonlinear  system  described  as  follows: 

xc  =  Acxc  +  fc{xc,xs)  =  Acxc  +  f?(xc,xs)  +  f*(xc,xs)  +  ---,  (7.3) 
xs    =    Asx,  +  fs(xc,  xs)  =  A'xs  +  tf(xc,xs)  +  f*(xc,  *,)  +  •••.  (7.4) 

Suppose  that  after  applying  an  appropriate  transformation  of  form 

xe  =  xc  +  t>i(*)  +  •••  +  «?(*)  (7-5) 
xs   =   x.  +  t»|(x)  +  •••+ «*(*),  (7-6) 

the  dynamics  for  (xc,  x3)  takes  the  form 

ic  =  Acxc  +  5c2(xc)  +  5c3(xc)  +  ---  +  ^(xc)  +  0(|^|'t+1)  (7.7) 
I,    =    A%  +  ^(ic,is)  +  5^(xc,xs)  +  .--  +  ^(ic,xs)  +  0(|x|,:+1),  (7.8) 

where  equation  (7.7)  has  the  exact  same  dynamics  as  the  critical  dynamics  of  a  normal  form 
(up  to  degree  k)  for  (7.3)  and  (7.4),  but  equation  (7.8)  does  not.  For  a  CASE  1  critical 
nonlinear  system,  the  dynamics  for  xc  is 

xc=akxk  +  0(\x\k+1)  (7.9) 

where  a*  is  the  first  nonzero  real  number  in  the  sequence  {02,03,04,  •  •  •}. 
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Similarly,  for  a  CASE  2  critical  nonlinear  system,  the  dynamics  for  xc  takes  the  form 

xcl    =s    -xc2  +  (x2cl  +  x2c2)(m3xci  -  n3xc2)  +  ■  ■  ■ 

+(*cl  +^c2)J(m(2j+i)^ci  -  n(2j+1)ic2)  +  0(\xc,xc2\2l+2)  (7.10) 
x~c2    =      uSd  +  {x2cl  +  x2c2)(m3xc2  +  n3xci )  +  •  •  • 

+(x2cl  +  x2c2)j(m{2j+1)xc2  +  nW+i)Xcl)  +  0(|xc,xs|(2j+2)),  (7.11) 

where  A;  =  2j  +  1  is  odd,  and  is  the  first  nonzero  real  number  in  the  sequence 
{m3,  m5,  m7,  ■■■}. 

We  will  prove  the  following  theorem  by  the  explicit  construction  of  Lyapunov  function 
candidates. 

Theorem  28  For  any  CASE  1  system  described  by  (7.9),  (7.8),  it  is  LAS  when  <  0  and 
k  is  odd,  and  unstable  otherwise.  For  any  CASE  2  system  described  by  (7.10),  (7.11)  and 
(7.8),  it  is  LAS  when  m2j+i  <  0,  and  unstable  otherwise.  OOO 

Before  proving  this  theorem,  we  need  the  following  lemmas. 

Lemma  8  Let  a  homogeneous  polynomial  function  Uk(x)  :  Rn  —>  R  be  positive  definite, 
where  k  >  0  is  an  even  integer.  Then  any  analytic  function  U(x)  :  Rn  — ►  R  defined  by 

Il(x)  =  Ilk(x)  +  0(\x\k+1),  (7.12) 


is  also  locally  positive  definite. 


000 
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Proof:  Let  fi  be  a  neighborhood  of  the  origin.  Let  0  /  i  6  0,  and  /  g  R  with  0  <  |/|  <  1. 
By  Taylor's  theorem  there  exists  a  real  parameter  0  <  £  <  1  such  that 

Consider  the  function  T(t)  =  ^^]j^\y=tx(x)k ■  The  function  T(t)  is  locally  continuous 
since  II(x)  is  Cw.  Since  T(0)  =  ^DkU(0)(x)k  =  Ilk{x)  >  0,  by  continuity  there  exist  a  real 
number  e  >  0  such  that  for  all  \t\  <  e,  T(t)  >  0.  This  means  that  for  all  \t\  <  min{l,£:}, 
^DkU(tx)(x)k  >  0.  Since  0  <  f  <  1,  we  have  ^DkU(^tx)(x)k  >  0.  Because  k  is  an  even 
number,  ^DkIl(^tx){x)k  =  U(tx)  >  0  whenever  0  and  x  £  0.  Thus  U(x)  is  locally 
positive  definite.  OOO 

Lemma  9  Consider  a  linear  system  x  =  Ax.  If  it  is  asymptotically  stable,  andV(x)  =  xTPx 
is  a  Lyapunov  function,  then  for  any  positive  integer  k,  W(x)  =  (xTPx)k  is  a  Lyapunov 
function  for  any  analytic  nonlinear  system  x  =  Ax  +  0(\x\2). 

Proof:  The  positive  definiteness  of  W(x)  is  evident.  Differentiating  W(x)  respect  to  time 
yields 

^  =  k{xTPxf-*^  =  k(xTPx)k-*  (xT(ATP  +  PA)x  +  0(\x\*)) 
at  at  '  ' 

Since  (xTPx)k~'1  is  positive  definite  and  is  locally  negative  definite,  dWW  \5  locally 
negative  definite.  OOO 


We  are  now  ready  to  prove  the  main  theorem  in  this  chapter. 
Proof  of  Theorem  28:  Since  a  CASE  2  system  is  equivalent  to  a  CASE  1  system  after  using 
polar  coordinates,  theorem  (28)  will  be  only  proved  for  CASE  1  systems. 
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•  ak  <  0  and  k  is  odd:  Choose  V(  xc,  x3)  =  ^x^  +  -fc  (XJP  x«)  ^  +  ^  which  is  positive 
definite.  Then 

£SSfi*sl  =  xc  (a*i*  +  0(|f  I**1))  +  (^x.)M/J  (fftf.  +  0<|.|»)) 

=  M*"  +  (xrPxs)(fc"1)/2  (fftf.)  +  0(|*|*+2). 

where  Q  =  (AsrP  +  PA4).  The  local  negative  definiteness  of  dV^x')  follows  since 
Tlk+1(xc,xs)  =  akxkc+1  +  (xJPx,)  *  1)/2  (xjQxs}  is  negative  definite. 

•  a*  >  0  and  k  is  odd:  Choose  F(  xc,xs)  =  -|*c  +  lE+T  (z^P^s)^1^2  which  is  locally 
indefinite. 


is  locally  negative  definite. 
•  ak  ^  0  and  k  is  even:  Choose  V{  xe,xt)  =  akxc  +  \  (xJPx^  1  which  is  indefinite. 

=  -41}  +  (xfpx,)"2-  (fr«i.)  +  o<|i  |*+') 


is  locally  negative  definite. 

These  observations  prove  Theorem  28  for  CASE  1  systems. 

For  CASE  2  systems,  the  local  asymptotic  stability  or  instability  can  be  proved  by  using 
V(  xc,xs)  =  ±\{xl,  +  xl2)  +  fc  (xJPxs)ik+1)/\  OOO 

Lyapunov  functions  have  been  constructed  for  CASE  1  or  CASE  2  nonlinear  systems 
in  partial  normal  forms.  One  of  the  properties  of  these  Lyapunov  functions  is  that  their 
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time  derivatives  along  the  system  dynamics  are  negative  definite  homogeneous  polynomial 
functions  plus  terms  of  higher  degree.  Lyapunov  functions  in  the  original  coordinates 
(xc,  xs)  can  be  obtained  by  just  substituting  the  transformation  x  =  v(x)  which  was  used 
to  obtain  the  partial  normal  forms,  i.e., 

V{xc,  Cg)  —  V(xc,  Xs)  \x=zv(x)  ■ 

In  the  next  section,  we  will  propose  a  method  for  simplification  of  Lyapunov  functions  in 
the  original  coordinates. 

7.3    Simplification  of  Lvapunov  Functions 

In  this  section  we  address  how  to  simplify  the  Lyapunov  functions  constructed  using 
the  algorithm  proposed  in  Section  7.2. 

For  an  LAS  CASE  1  nonlinear  system,  a  Lyapunov  function  in  partial  normal  form 
coordinates  (xc,xs)  constructed  in  previous  section  is  given  by 


F(xe,i,)  =  |x2  +  -lT(ffP^)(i+1) 


where  a2k+i  is  the  first  nonzero  real  number  in  the  sequence  {a2,  03,  a4,  •  •  •}. 

Similarly,  for  an  LAS  CASE  2  nonlinear  system,  a  Lyapunov  function  in  partial  normal 
form  coordinates  (xc,xs)  constructed  in  previous  section  is  given  by 

1  /-rD- 


V(i c,  x.)  =  i(x2,  +  ij2)  +         (xJP  xs) 


where  m2k+\  is  the  first  nonzero  real  number  in  the  sequence  {1713,  05,  07,  •  •  •}. 
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A  Lyapunov  function  in  the  original  coordinates  (xc,  xs)  is  obtained  by  just  substituting 
the  forward  transformation 


x  =  v(x)  =  x  +  v2(x)  +  •■•  +  v2k+1(x), 

into  V(x)  to  give 



V(xc,  xs)  =  y(ic,i,)|x=v(x) 
=    V\xc,  x.)  +  ■  ■  ■  +  V*k+2(xc,  xt)  +  •  •  •  +  V(2k+W+1\xc,  xs),  (7.14) 

which  is  quite  complicated,  and  includes  terms  of  degree  2  to  (2k  +  l)(k  + 1).  Consider  now 
the  time  derivative  of  (7.14)  along  the  dynamics 

dvix^  =  ^£^^)  =  (n^(ia,,.)+o(|f|^))^ 

=   n2k+2(xeix,)  +  0(\x\2k+%  (7.15) 

for  some  homogeneous  polynomial  function  of  degree  2k  +  2:  Yl2k+2(xc,  xs).  Notice  that 
term  V2k+3(xc,  &,)  +  •••  +  V^2k+1^k+1\xc,  xs)  in  (7.14)  only  contributes  terms  of  degree 
2k  +  3  or  higher  in  (7.15).  Let  function  Vnew(xc,  xs)  be  obtained  by  truncating  the  terms 
of  degree  2k  +  3  or  higher  in  (7.14),  i.e., 

Vnew(xc,  x.)  =  V2(xc,  xs)  +  V3(xc,  x.)  +  ■  ■  ■  +  V2k+2(xc,  x,).  (7.16) 
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The  time  derivative  of  Vnew(xc,  xs)  will  also  take  the  form: 


The  following  theorem  is  obtained  due  to  this  observation. 

Theorem  29  Assume  that  a  CASE  1  or  CASE  2  critical  nonlinear  system  x  =  f(x)  is 
LAS,  and  a2k+i  <  0  (or  m2k+i  <  0)  is  the  first  nonzero  real  number  in  the  sequence 
{a2,  a3,  a4,  •  •  •}  (or  {m3,  m5,  m7,  •■■})■  IfV(xc,  xs)  constructed  as  in  (7.14)  is  a  Lyapunov 
function  for  x  =  f(x),  then  Vnew(xc,  xs)  =  V2(xc,  xs)  +  V3(xc,  !,)  +  •••+  V2k+2(xc,  xs) 
is  a  Lyapunov  function.  OOO 

Proof:  The  local  negative  definiteness  of  dVnew(xc,  xs)/dt  is  already  proved.  We  need 
to  prove  that  Vnew(xc,  xs)  is  locally  positive  definite.  We  will  prove  it  by  contradiction. 
Remember  that  locally,  a  function  can  only  be  indefinite,  positive  semi-definite,  positive 
definite,  negative  semi-definite,  or  negative  definite. 

•  Assume  that  V"ne^(ic,  xs)  is  locally  indefinite  (or  negative  definite,  or  negative  semi- 
definite).  Therefore,  there  exists  a  region  arbitrarily  close  to  the  origin  in  which 
both  dVnew(xc,  x,)/dt  and  Vnew(xc,  xs)  are  negative.  This  implies  that  the  system 
is  unstable  using  an  instability  theorem  given  in  [37].  That  is  a  contradiction. 

•  Assume  that  Vnew(xc,  xs)  is  positive  semi-definite  and  dVnew(xc,  x„)/dt  negative  in 
an  open  neighborhood  of  the  origin  (xc  =  0,  xs  =  0).  Then  there  exists  a  point  x°  € 
Q  such  that  Vnew(x®,  x°)  =  0.  By  assumption,  there  is  an  open  neighborhood  A  C  fi 
of  x°  such  that  for  any  y  £  A  C  ft,  Vnew(yc,  ys)  >  0.  However  dVnew(yc,  y,)/dt  <  0 
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implies  that  there  exists  an  point  z  £  A  C  SI  such  that  Vnew(zc,  zs)  <  0.  This  is  alsoy 
a  contradiction. 

In  summary,  Vnew(xc,  x.)  =  V2(xc,  xs)  +  V3(xc,  xs)  +  ■■■  +  V2k+2{xc,  xt)  must  be 
locally  positive  definite,  and  hence  a  Lyapunov  function.  OOO 

7.4    Computation  of  Partial  Normal  Forms 

In  this  section,  we  address  the  issues  related  to  the  computation  of  partial  normal 
forms  and  the  associated  transformations.  Let  the  original  system,  the  transformed  system, 
transformations  and  the  PDE  be  given  respectively  as  follows: 


/  \ 

i. 


\  Xs  ) 
I  \ 

/  \ 

Xc 


(  \ 

Vc 


(^ 

da 
du, 

V  9a  / 


fc(x,a 
f3(x,a 

9c(y,a 

gs(y,a 
uc(y,a 

us(y,a 

vc(x,  a 
vs(x,a 

Uc(u,a 
Us(u,a 


\  I 

oo 


m=0 


f?+\x) 

\  /r+1(*) ; 
/  \ 

5cm+1(y) 


ml ' 


m=0 


g?+\y) 


ml ' 


m=0 


^  <+1(y)  ) 


m 


,  u(y,0)  =  y, 


( 


m=0 


\ 


\  I 

CO 


m=0 


J 


U?+2{u) 
U?+\u) 


m\ 


;,  V(X,Q)  =  X, 


m: 


-,  u(y,0)  =  y. 


(7.18) 


(7.19) 


(7.20) 


(7.21) 


(7.22) 


The  terms  of  degree  m,  gm,  fm,  and  Um  are  related  by  the  following  equation,  where 
Tm  only  depends  on  /\  •  •  • ,  fm  and  U2,  ■  ■  ■ ,  Um-\ 


gm(x)    =  Tm{x)-LJVm{x) 
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=  rm(x)-( 


dUm(x) 


dx 


Ax  -  AW 


'(*)) 


17(0!) 


{ /  \  / 

'  aiv  at/,."* 


-  < 


9xc  9xa 


Acxc 


\ 


) 


acu: 


ASU? 


I J 


dU, 


It  is  clear  that  the  normalization  of  g™(x)  depends  on  the  choice  of  U™(x).  Therefore  we 
can  simply  set  U™(x)  =  0  to  simplify  the  computation  of  partial  normal  forms.  Using  this 
construction,  Theorem  12  and  Theorem  13  in  Chapter  3  can  be  reduced  to  the  following 
forms. 

Theorem  30  Let  f(x,a),  g(y,a)  and  U(u,a)  be  expanded  as  above,  and  take  fi  (or  gi,  or 
Ui)  to  be  f,+1  (or  g,+i ,  or  Ut+2),  respectively.  The  sequence  {ScA^Kii  may  ^6  computed 
as  follows: 

( 

9c,i(v) 

9s,i{y) 


\ 


f&tv) 


I 


\ 


(rn) 

where  the  sequence  f\_^  is  defined  recursively  by 


I  ,<»>  N 

Jc,t 

AO) 


/  \ 

fc,i 


and  for  m  =  1, 2, 3,  •  •  • ,  i, 


f(m)  (x) 


r*(n?_1)  (x) 

/"(m_1)  (x) 
^  Js^—m+iy*)  J 
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i  -  m 


dxc 


'  IT     /-'v      Qt/c.jtx)  /   x      9t/c,,(i)  f(m-l)     ,  v 

^c.jW  gXc     JC),-_m_jVx  J  9ia  Js,i-m-j\x) 


\  dxc  uc,3\J-J 


Theorem  31  Assume  u(y.a)  is  the  solution  of  the  PDE,  and 


000 


I 


p{x,a)  = 


pc(x,a) 
ps(x,a) 


I 


E 

m=0 


Pc,m(x) 
Ps,m{x) 


\ 


TO! 


(7.23) 


Let 


The  sequence  {<?^,,  <?J,}  may  be  computed  as  follows: 


I  \ 

qdy,a) 


9s(y,a) 


( 


=  p(u(y,a),a)= 


771  =  0 


Qc,m(y) 
Qs,m(y) 


\ 


(7.24) 


9.(2/)  = 


/  \ 

9cAy) 
Qs,i{y) 


p%(y) 
p%(y) 


where  the  sequence  p^™)m  is  defined  recursively  by, 


(o) 


(m)  /  \ 

Pi-m(X)  = 


1  (o)\ 

Pc/i 

1  P(0)- 
\  Ps>>  ) 


(  \ 

Pe,i 


=  Pi  = 


Km) 


{P(™i-m(x)  ) 


P{™i-m+l(x) 


\ 


+  £ 

J=0 


t  —  m 

i    3  j 


Ucj(x),  to  =  1,2, 


1  9P^iW  N 


9xc 
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000 

The  algorithms  for  deriving  the  transformations  um,  vm,  can  be  easily  written  out 
accordingly.  With  a  little  bit  more  effort,  we  can  easily  verify  that  u™  =  v™  =  0  if  £/cm  =  0 
for  m  =  2,3, The  first  few  elements  in  the  sequence  {u2c,  u*,  u*,  ■  ■  •}  and  {v],  vf,  v*, 
•  •  •}  can  be  written  out  as  follows: 

»?(»)  =  vKv) 

«Sd)  =  uH,)+^MuH,) 


and 


»c3(y)  =  -^^2(y)-^) 

^  =  _^{y)_2^{y)_uHy) 


Remark  10  The  advantage  of  the  computation  of  partial  normal  forms  and  the  associated 
transformations  can  be  observed  from  the  discussions  above.  Assume  x,  and  xc  have  the 
same  dimension.  The  amount  of  time  necessary  to  compute  partial  normal  forms  and  the 


109 


associated  transformations  is  about  half  of  the  time  necessary  for  computing  normal  forms 
and  the  associated  transformations  for  all  states.  Actually,  the  dimension  of  xc  is  much 


7.5  Examples 

In  this  section,  examples  are  provided  to  illustrate  the  use  of  the  algorithm  proposed  in 
this  chapter. 

Example  21  (Reconsideration  of  Example  23) 
Consider  again  the  system 


lower  than  the  dimension  of  xs  in  many  applications. 


OOO 


x 


x3  +  x2y2, 


(7.25) 


if 


y  +  x2  +  y3. 


(7.26) 


A  Lyapunov  function  is  given  by 


with 


dV(x,y) 


=  -(x4  +  y4)  +  0(\x,y\5). 


dt 


OOO 


Example  25  (The  CASE  1  system  in  Example  16  in  Chapter  5) 


Consider  the  system 


x  = 


-xy  +  xz2  +  xyz, 


y    =    -y  +  x2  +  xz2, 


z   —    —z  +  xy. 


A  partial  normal  form  (up  to  degree  3)  is 


y   =    -y  +  x2  +  xz2  +  2x2y, 


I    =    -z  +  xy  +  xy2. 


The  associated  backward  transformation  is 


xyz     xz2  xy2 
x    =    x  +  xy-  —  -  —  +  —, 


y  =  y, 


z    =  z. 


The  associated  forward  transformation  is 


xy2     xyz  xz2 
x    =    x-xy+  —  +  —  +  —, 


y  =  y, 


z   =  z. 


The  following  Lyapunov  function 


V(x,  y,  z)  =  x2  -  2  x2y  +  z2x2  +  2  x2y2  +  yzx2  +  y4  +  z4 
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Ill 


is  derived  by  truncating  terms  of  degree  of  5  or  higher  using  V(x,  y,z)  =  x2  +  y4  +  z4 .  000 
Example  26  (The  CASE  2  system  in  Example  17  in  Chapter  5) 
Consider  the  system 


x    =    —y  +  xz2  +  xyz  +  x2-\-z2, 


y   =    x  +  z2  +  i2  +  xz2, 


z    =    —z  +  xy  +  z2. 


A  partial  normal  form  (up  to  degree  3)  is 


.     x3     7x2y     xy2     7  y3 

x    =  —v  

y     4        12        4  12 

_    (3a -7y)  (x2  +  f) 

=  -y  5  

7  x3     x2y     7  xy2  y3 

y  =  x  - 


=    x  — 


12       4        12  4 
(7i  +  Sy)  (x2  +  y2) 
12 


x2y     Axy2     2y3     3yz2     x3  xz2 
z    =    -z  +  xy  +  z+  —  +  -  ___-  +  ___ 


The  associated  backward  transformation  (up  to  degree  3)  is 


£*_     2xy     2f_     Zz2     7 x3     x2y  12x2z 
X   ~   X  ~  3  +  ~3        3__5-36~2_  25 

Zlxyz     xz2     11  y3     12  fz     yz2     2z3  2xy2 
50~ +  10  ~   12   +    25    ~  20  ~   5   +    9  ' 
x2     2xy     2y2     z2     x3     5x2y     7  x2z 
y  -    y  +  Y  +  ~y~  +  T"~"5"+12~   18  +  50 
17  xyz     xz2     13  f     7y2z  2xy2 
50~ +  _20"~~36        50~+    3  ' 

z   =  z. 
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The  associated  forward  transformation  is 


x2     2xy     2y2     Sz2     23  x3     x2y  I2x2z 

x     -    x  4  -  4  — I  1   -I  

3       3        3        5        36        6  25 

2xy2     Zlxyz     xz2     by3     12  y2z  yz2  2z3 

+    9    +    50    +  ~6        12        25  12 +  T~1 

x2     2xy     2y^            x^_     25x2y  lx2z 

V   ~    V~  Z~   Z    ~   3   +y~12+    18  50 

2        17t,i,       7-r,2       OQ  „3       7„2,  O 


2  xy      Ylxyz     1  xz1     29  y3     7y'z  2yz 


+  "^  +  ^T  7^  +  ^7T  + 


3  50  12        36        50         3  ' 


z   —  z. 


The  following  Lyapunov  function 


a      ,     2x3        »      6xz2     4  y3     2yz2  3x4 
V(x,y,z)   =    x2  +  y2  +  —  -2x2</+  —  3   +    5    +  ~ 

24  zx3  yx*        2  2     24zx2y     3z2x2  xy3 

+    25  2  y         25  5  2 

7*xy2  12z2xy     4z3x     5y4     7  zy3     4  z2t/2  7z4 
25  5      +    5+   2   +    25  ~5+5 


is  derived  using  V(x,y,z)  =  x2  +  f  +  z4.  000 

7.6  Summary 

In  this  chapter,  an  alternative  method  is  proposed  for  the  construction  of  Lyapunov 
functions  for  CASE  1  and  CASE  2  critical  nonlinear  systems.  This  method  requires  deriving 
only  a  partial  normal  form  up  to  an  appropriate  degree  so  that  the  computation  is  reduced 
significantly.  After  such  a  partial  normal  form  is  obtained,  a  Lyapunov  function  for  the 
system  in  the  partial  normal  form  coordinates  is  constructed  easily.  A  Lyapunov  function 
for  the  system  in  the  original  coordinates  is  obtained  by  just  substituting  of  the  associated 
forward  transformation.  This  Lyapunov  function  can  be  further  simplified  by  truncating 
some  of  the  terms  of  higher  degrees. 


CHAPTER  8 
APPLICATIONS 

Two  different  algorithms  for  the  construction  of  Lyapunov  functions  for  CASE  1  and 
CASE  2  critical  nonlinear  systems  in  were  developed  in  Chapters  5  through  7.  In  Section 
8.1  of  this  chapter,  nonlinear  systems  with  multiple  critical  eigenvalues  are  studied  using 
normal  form  techniques.  This  dissertation  focuses  on  continuous-time  nonlinear  systems. 
Normal  form  techniques  are  also  applicable  to  discrete-time  nonlinear  systems  with  only 
slight  modifications.  An  introduction  to  normal  forms  for  discrete-time  nonlinear  systems 
is  given  in  Section  8.2.  Finally,  a  brief  summary  of  this  chapter  is  given  in  8.3. 

8.1    Nonlinear  Systems  with  Multiple  Critical  Eigenvalues 

In  this  section  we  study  critical  nonlinear  systems  with  more  than  one  zero  eigenvalue  or 
more  than  one  pair  of  purely  imaginary  eigenvalues  using  normal  form  techniques.  Sufficient 
conditions  for  the  LASY  of  such  nonlinear  systems  will  be  derived.  Let  us  look  at  an  example 
first. 

Example  27  Consider  a  critical  nonlinear  system  as  follows: 


X    =    -y  -  y2  +  2xy  -  4x3  -  ^x2y  -  2xy2  -  ^y3 

-4  x4  +  -x  V  +  xy3  -™y*-4x5-  —y5 
2  2  12 

+|x42/  +  4  x3y2  +  Yx*y3  + 7  xv4 

y   =    x  +  x2  -  2  xy  +  |x3  -  6  x2y  +  ^  xy2  -  4  y3 


(8.1) 
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2  = 


+  *LX4  _  17  x3y  +  H  xy  _  4  4  +  323  5 
2  *     2  12 

-4  y5  -  39  x4j/  +  -^-z V  +  4  x2j/3  -  |  x j/4 

9        2  2      *^  2  2  2        ^  y 

y  —  x  —  xyi  —  +  xz-yz  +  xy  +  — 

lli4     y2z2     5y3z     2x2yz     19  xy3 
fi  3  3  3 


2 

11  u4     29x3u     2xy2z     x2z2  5x3z 

— —  +   +  — —  +  +  

2  3  3  6  3 


+ 


119x5 

x2yz2 

+ 

25  xV 

8 

12 

8 

119  y5 

43  xy3 

53  x3yz 

8 

6 

6 

7y3z2 
24 

2  2 

xy  z 
12 

+ 

7xV 
24 

12 


r3,,2 


71  j4; 
12 


(8.2) 


(8.3) 


The  Jacobian  matrix  of  this  nonlinear  system  at  the  origin  is 


/  \ 
0-10 


1    0  0 


i  0    0    0  j 


.  This  ma- 


trix possesses  a  zero  eigenvalue  and  a  pair  of  purely  imaginary  eigenvalues.  Without  any 
simplification,  this  system  looks  very  complicated.  However,  if  the  transformation 


x  = 


V  = 


x  -  x"  +  -x% 


1 


y  y  1  ^y  ^y 

,     x  y     xyi     3r  xy  z 


2,2 


2x2y2     bxyz2     xlyz     xyz6     z      7  xy2z 

3  6  3~~ +  ~6       ~6~  +  24 

x4y  +  7x3y2     x3yz  |  7x2y3  |  lx2y2z 


24 

8 


7 xlyz 
24~ 


12 
xy32 
12 


24 


12 


(8.4) 
(8.5) 


(8-6) 
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is  applied  to  (8.1)-(8.3),  we  obtain  a  new  system  as  follows: 

x    =    -y-  (x2  +  y2)(3y  +  4x)  +  0(\-\4)  (8.7) 

g   =         (x2  +  y2)(Zx-4y)  +  0(\-\4)  (8.8) 

I   =    -z5  +  0(\-\6).  (8.9) 


If  we  take 

V(x,y,z)  =  (x2  +  y2)2  +  z2  (8.10) 

as  a  Lyapunov  function  candidate  for  the  system  (8.7)-(8.9),  it  is  straightforward  to  verify 
that 

=  _16  (  x2  +  f  )3  -  2  t  +  0(|  •  |7).  (8.11) 

at 

The  LASY  of  the  system  (8.7)-(8.9)  (or  equivalently,  that  of  (8.1)-(8.3))  follows  since  (8.11) 
is  locally  negative  definite.  A  Lyapunov  function  in  the  original  coordinates  (x,  y,  z)  can  be 
derived  by  just  substituting  equations  (8.4)-(8.6)  into  (8.10).  OOO 
From  this  example,  we  can  make  the  following  observations.  As  before,  the  critical 
nonlinear  systems  considered  here  are  supposed  to  have  following  structure: 

xc  =  Acxc  +  /c2(xc,  !,)  +  ••■  (8.12) 
xa    =   A3xs  +  f2(xc,xs)  +  ---  (8.13) 

where  all  of  the  eigenvalues  of  A*  are  stable,  while  all  of  the  eigenvalues  of  Ac  are  critical. 
By  choosing  an  appropriate  transformation  of  the  the  form: 


xc  = 


xc  +  v2(xs,  xe)  +  •  •  •  +  v*(xa,  xc) 
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xs    =    x3  +  v]{xs,xc)  +  ■  ■  •  +  w*(xs,ic) 

the  dynamics  for  (xc,  xs)  satisfy 

le  =  Acxc+g2c(xc)  +  ---  +  gkc(xc)  +  0(\xc,xs\k+1)  (8.14) 
Sa    =    Asis  +  <7,2(xc,xs)  +  ---  +  ^(^,is)  +  0(|ic,ia|fc+1).  (8.15) 

Theorem  32  Assume  that  the  subsystem  (8.14)  *s  LAS  without  the  necessity  of  considering 
the  higher  order  terms,  and  furthermore  a  Lyapunov  function  V(xc)  exists  such  that 

9  («.)  =  +  9l{*c)  +  ■■■  +  9kc(xc))  (8.16) 

is  a  homogeneous  polynomial  of  degree  2r.  Let  P  =  PT  >  0  s«c/i  </iat  V(xs)  =  xj  Pxs  is 
a  Lyapunov  function  for  the  linear  system  is=  Asxs.  Under  these  hypothesis,  the  function 
V(xc,xs)  =  (xj  Pxsy  +  V(xc)  is  a  Lyapunov  function  for  the  system  (8.14)-(8.15).  OOO 

In  [23],  Fu  has  also  done  some  research  on  the  construction  of  Lyapunov  functions  for 
critical  nonlinear  systems  with  multiple  eigenvalues  when  A°  is  diagonalizable.  That  work 
is  similar  to  what  is  presented  below  in  Corollary  1  except  for  the  fact  that  it  can  only 
be  applied  to  critical  systems  for  which  the  LAS  can  be  determined  by  the  first  through 
third  order  terms.  After  those  systems  are  transformed  into  some  sort  of  normal  forms, 
the  construction  of  Lyapunov  functions  for  these  classes  of  critical  nonlinear  systems  can 
be  easily  understand.  In  fact,  the  algorithm  presented  in  Chapters  7  for  the  construction 
of  Lyapunov  functions  for  CASE  1  and  CASE  2  critical  nonlinear  systems  is  only  a  special 
case  of  Theorem  32. 
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Corollary  1  Consider  a  critical  nonlinear  system  of  the  form  (8.12)  and  (8.13)  with  (without 

of  loss  of  generality)  Ac  =  diag(A\,  A2,  •  •  • ,  A„c),  where  for  j  =  1,2,  ■•■,nc,  Aj  is  either 

(  \ 


a  scalar  zero  or  Aj  = 


0  -ui 


\ 


with  Uj  €  R+  ■  If  xcj  (j  =  1,2,  •••,nc)  is  the  state 

u3  0 

cooresponding  to  Aj,  and  the  critical  dynamics  xcj  in  a  (partial)  approximate  normal  form 


of  the  form  (8.14)  and  (8.15)  satisfy 


2Jfc,+l 


lci=  £  artj  +  0(\x\2k>+2) 


i-2 


(8.17) 


if  Aj  =  0  where  aJi2  =  •  •  •  =  a,^  =  0  but  aj,2fc,+i  <  0, 


or 


lcjil    =    -UjXejt2  +  Y^imjii+iXcji  +  nji2i+lxCh2)  +  0{\x\2k>+2)  (8.18) 
1=1 

icj,2    =   UjXcji  +  ^(rnjii+iXcja-nj,2i+iXcjA)  +  0{\x\2ki+2)  (8.19) 


:=1 


0  -Ui 


\ 


if  Aj  =  where  mj,3  =  •  •  •  =  mj)2fc-i  =  0  but  mj^k}+i  <  0,  then  the  nonlinear 

Uj  0 

system  (8.14)  and  (8.15)  (or  equivalently  (8.12)  and  (8.13))  is  LAS. 


8.2    Normal  Forms  For  Discrete- Time  Nonlinear  Systems 

In  this  section,  an  introduction  to  normal  form  techniques  for  discrete-time  nonlinear 
systems  is  provided. 

Consider  a  discrete-time  system  represented  as  follows, 


x  •  ►  f(x)  =  Ax  +  f2(x)  +  ■■■  +  fk(y)  + 


(8.20) 
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Here  y->  represents  iterative  transition.  Let 

x  =  y  +  Pk(y) 

be  a  near  identity  change  of  coordinates.  The  dynamics  for  y  will  take  the  form 

y~Ay  +  g\y)  +  ■■■  +  fk(y)  -  {Pk(Ay)  -  APk(y)}  +  gk+\y)  +  ■■■  (8.21) 

where  gi  =  /•'  for  i  =  2,3,- 1,  gk(y)  =  fk(y)  -  {Pk{Ay)  -  APk(y)},  and  gi  ?  f  in 
general  for  i  >  k  +  1.  Therefore,  the  linear  operator  L\  :  H(n,n,k)  — ►  H(n,n,k)  can  be 
defined  by 

LkAPk(y)  =  Pk(Ay)-APk(y).  (8.22) 

For  discrete  time  systems,  a  eigenvalue  A  is  called  critical  if  and  only  if  its  magnitude 
11*11  =  1. 

Definition  8  [4]  Let  a(A)  =  {Aj,  A2,  An}  be  the  spectrum  of  the  n  x  n  matrix  A.  Let 
a  =  [c*i,  •  •  -  ,a„]  be  a  set  of  nonnegative  integers.  A  monomial  xaej,  with  5Z"=i  =  ^  >  2 
is  called  a  resonant  monomial  if 

Aa  dLf\^...\%»  =  X3  (8.23) 

for  some  integers  k  and  j.  This  condition  is  called  the  resonant  condition  for  the  matrix  A 
in  discrete-time  domain.  000 

One  of  the  fundamental  theorems  about  normal  forms  for  discrete-time  systems  is  given 
as  follows. 
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Theorem  33  [4]  If  A  is  an  upper  triangular  with  (ordered)  diagonal  elements  {Ai,  A2,  •  •  • ,  An}, 
then  a  normal  form  map  can  be  chosen  such  that  its  nonlinear  part  contains  only  resonant 
monomials.  000 

Without  any  difficulty,  one  can  prove  the  following  theorem. 

Theorem  31  Consider  a  discrete-time  critical  nonlinear  system: 


(  \ 

I 

\  Xs ) 

Acxc  +  f?(xc,  x.)  +  ff(xe,  xs)  + 
A°xs  +  f?(xc,xs)  +  f?(xc,xs)  + 


(8.24) 


where  all  of  the  eigenvalues  of  Ac  are  critical  while  all  of  the  eigenvalues  of  A3  are  stable. 
Any  normal  form  for  (8.24)  UP  t°  degree  k  has  the  following  structure: 

\ 


(  \ 

1 

Xc 

V  is  ) 

A<xc  +  Pc{xc)  +  ■■■+  +j»(*e)  +  0(xc,  */+*), 


where  f's(xc,  xt  =  0)  =  0  for  i  =  2, 3,  •  •  •,  k. 


(8.25) 


8.3  Summary 

In  this  chapter  we  first  discussed  how  to  construct  Lyapunov  functions  for  LAS  nonlinear 
systems  with  multiple  critical  eigenvalues.  Following  this,  an  introduction  to  the  normal 
form  techniques  for  discrete-time  critical  nonlinear  systems  was  presented. 


CHAPTER  9 
SUMMARY  OF  THE  DISSERTATION 

Nonlinear  dynamical  models  arise  from  both  theoretical  research  and  real  applications, 
including  robotics,  aircraft  control,  adaptive  control,  chemical  processes,  numerical  algo- 
rithms, and  power  systems.  The  most  important  analysis  for  nonlinear  systems  is  the 
determination  of  stability.  As  far  as  stability  is  concerned,  only  critical  nonlinear  systems 
need  further  investigation.  For  critical  nonlinear  systems,  stability  analysis  is,  up  to  now, 
predominantly  performed  using  center  manifold  and  bifurcation  theories  in  combination 
with  Lyapunov  techniques.  When  a  dynamical  system  is  asymptotically  stable,  a  Lya- 
punov  function  must  exist.  Many  system  properties,  such  as  stability  region,  convergence 
rate,  and  robustness,  are  associated  with  and  can  be  estimated  using  Lyapunov  functions. 
The  search  on  the  determination  of  stability  properties  of  critical  nonlinear  systems  or  the 
construction  of  the  associated  Lyapunov  functions  is  far  from  complete. 

In  this  dissertation,  critical  nonlinear  systems  are  studied  using  normal  form  theory. 
Once  a  critical  nonlinear  system  is  in  a  normal  form,  the  dynamics  of  the  critical  states  is 
decoupled  locally  from  those  of  the  stable  states.  Similar  to  results  derived  using  center 
manifold  theory,  the  stability  of  the  entire  nonlinear  system  is  determined  by  that  of  the 
critical  dynamics.  Thus  the  normal  form  method  is  also  a  model  reduction  technique.  The 
relationship  between  two  model  reduction  techniques,  i.e.,  normal  form  and  center  manifold 
methods,  is  explored.  Compared  to  the  center  manifold  method,  the  normal  form  method 
provides  not  only  a  model  reduction  technique,  but  also  a  direction  for  the  construction 
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of  Lyapunov  functions.  If  the  critical  dynamics  is  locally  asymptotically  stable,  and  a 
Lyapunov  function  is  available,  then  Lyapunov  functions  for  the  entire  system  can  be  easily 
constructed  as  a  sum  of  two  indivadule  Lyapunov  functions  corresponding  to  the  critical 
and  stable  dynamics,  respectively. 

The  convergence  of  normal  forms  for  a  given  nonlinear  system  is  not  always  guaranteed. 
Construction  of  Lyapunov  functions  for  purely  critical  nonlinear  systems  is  only  known  in 
a  few  special  cases.  CASE  1  and  CASE  2  critical  nonlinear  systems  are  studied  intensively 
using  approximate  normal  forms.  The  explicit  stability  criteria  as  well  as  explicit  algorithms 
for  construction  of  Lyapunov  functions  are  proposed  once  such  a  system  is  put  into  an 
appropriate  (approximate)  normal  form.  However,  it  is  often  desirable  to  have  Lyapunov 
functions  for  the  original  systems.  By  further  studying  the  normal  form  transformations, 
it  is  shown  that  Lyapunov  functions  in  the  original  coordinates  can  be  derived  by  just 
substituting  an  appropriate  transformation. 

When  compared  with  Lyapunov  functions  constructed  by  other  researchers  for  some 
very  simple  critical  nonlinear  systems,  the  Lyapunov  functions  we  constructed  are  quite 
complicated  in  some  sense.  Methods  for  simplification  of  Lyapunov  functions  are  also 
developed  in  this  dissertation.  Using  the  results  in  Chapter  6,  the  Lyapunov  function  we 
derived  can  be  as  simple  as  those  constructed  by  other  researchers. 

The  derivation  of  normal  forms  and  the  associated  transformations  is  very  computation- 
ally intensive.  The  issue  of  reducing  the  computational  complexity  is  studied.  In  Chapter 
7,  an  alternative  algorithm  for  the  construction  of  Lyapunov  functions  is  proposed  wherein 
only  partial  (approximate)  normal  forms  are  necessary.  The  simplification  of  the  Lyapunov 
functions  is  also  addressed  in  Chapter  7.  A  more  general  of  class  of  critical  nonlinear  sys- 
tems is  examined  in  Chapter  8,  and  some  sufficient  conditions  for  determining  the  local 
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stability  and  for  constructing  Lyapunov  functions  are  obtained.  Finally,  the  use  of  normal 
form  techniques  for  discrete-time  nonlinear  systems  analysis  is  approached  in  Chapter  8. 
All  of  the  results  obtained  in  this  dissertation  for  continuous-time  critical  nonlinear  systems 
will  carry  over  to  discrete-time  with  some  modifications. 
In  conclusion,  the  following  has  be  accomplished. 

•  Equivalence  of  center  manifolds. 

•  Decoupled-like  structure  of  normal  forms  for  critical  nonlinear  systems. 

•  Relationship  between  center  manifold  and  normal  form  theories. 

•  Stability  analysis  of  critical  nonlinear  systems  using  normal  forms. 

•  A  general  direction  for  the  construction  of  Lyapunov  functions  for  critical  nonlinear 
systems 

•  Two  different  explicit  algorithms  for  the  construction  of  Lyapunov  functions  for  crit- 
ical nonlinear  systems  in  two  special  cases. 

•  Simplification  of  Lyapunov  functions. 

•  Simplification  of  the  computation  for  the  construction  of  Lyapunov  functions,  normal 
form  techniques. 

•  Introduction  to  the  use  of  normal  form  techniques  for  discrete-time  critical  nonlinear 
systems. 
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